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ABSTRACT 


— ~^k  linear,  unconstrained  perturbation  model  for  the  Tracked  Ram  Air 
Cushion  Vehicle  is  developed.  This  model  is  the  result  of  theoretical 
expressions  for  the  TRACV  which  have  been  verified  by  wind  tunnel  and 
towed  model  tests.  This  model  is  varied  to  allow  for  passively  suspended, 
two-degree -of- freedom  winglets  and  for  processor-controlled  actuators  on 
the  same  winglets.  Optimization  of  the  springs  and  dampers  in  the  passive 
suspension  is  performed  according  to  a performance  index  based  on  acceler- 
ation, winglet  gap  variation,  and  control  power.  Linear  optimal  control 
is  applied  to  the  active  suspension  to  determine  the  optimal  feedback 
gains  using  a similar  performance  index.  The  basic,  passively  suspended, 
and  actively  suspended  vehicles  are  analyzed  to  determine  root  mean  squared 
values  for  the  following:  l)  vertical  acceleration  in  the  foremost  and 
rearmost  seats  in  the  passenger  cabin,  2)  gap  variation  at  the  front  and 
rear  winglet  areas,  and  3)  control  deflection.  The  acceleration  spectral 
density  of  each  of  the  vehicle  types  is  compared  to  the  Urban  Tracked  Ram 
Cusion  Vehicle  standard.  The  active  control  system  is  analyzed  to  see  if 
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ABSTRACT 


__^»A  linear,  unconstrained  perturbation  model  for  the  Tracked  Ram  Air 
Cushion  Vehicle  is  developed.  This  model  is  the  result  of  theoretical 
expressions  for  the  TRACV  which  have  been  verified  by  wind  tunnel  and 
towed  model  tests.  This  model  is  varied  to  allow  for  passively  suspended, 
two-degree-of-freedom  winglets  and  for  processor-controlled  actuators  on 
the  same  winglets.  Optimization  of  the  springs  and  dampers  in  the  passive 
suspension  is  performed  according  to  a performance  index  based  on  acceler- 
ation, winglet  gap  variation,  and  control  power.  Linear  optimal  control 
is  applied  to  the  active  suspension  to  determine  the  optimal  feedback 
gains  up-'-"  a similar  performance  index.  The  basic,  passively  suspended, 
an<!  suspended  vehicles  are  analyzed  to  determine  root  mean  squared 

Vf  uiie  following:  l)  vertical  acceleration  in  the  foremost  and 

rearmost  seats  in  the  passenger  cabin,  2)  gap  variation  at  the  front  and 
rear  winglet  areas,  and  3)  control  deflection.  The  acceleration  spectral 
density  of  each  of  the  vehicle  types  is  compared  to  the  Urban  Tracked  Ram 
Cusion  Vehicle  standard.  The  active  control  system  is  analyzed  to  see  if 


a reduced  set  of  sensors  may  achieve  acceptable  ride  quality  based  on  the 
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CHAPTER  I 
INTRODUCTION 

Reliable  and  inexpensive  high-speed  intercity  transportation  is 
receiving  increasing  attention  as  the  cost  of  operating  automobiles  con- 
tinues to  soar.  Present  modes,  however,  face  time  and  convenience  problems 
which  the  car  does  not.  Aircraft  are  capable  of  cruising  at  great  speed, 
but  the  time  lost  on  the  ground  getting  into  and  out  of  airports  can  be 
thrice  the  enroute  time  on  short  and  medium  distance  runs.  In  addition, 
due  to  space  requirements , airports  usually  cannot  be  located  in  or  near 
the  center  of  cities.  Trains  do  not  face  the  same  requirements;  their 
stations  are  often  in  the  heart  of  the  business  district  of  town.  But 
while  trains  have  terminal  convenience,  they  often  have  as  little  as  one- 
tenth  the  cruise  speed  of  aircraft. 

The  most  convenient  form  of  mass  transit  would  combine  the  better 
h»lf  of  each  of  these,  giving  both  high  speed  and  terminal  convenience  to 
the  vehicle.  High  speed  ground  transportation  can  provide  this  combination 
once  the  dangers  to  public  safety  and  construction  problems  are  removed. 
Several  schemes  have  been  advanced  and  the  most  likely  fall  into  three  major 
categories:  high-speed  rail,  magnetic  levitation  (mag-lev),  and  air  cushion 
vehicles.  Each  of  the  classes  has  varying  subtypes  with  advantages  and 
drawbacks . 

The  high-speed  rail  concept  utilizes  an  improved  rail  with  welded 
joints  and  modified  existing  equipment  and  technology.  On  the  other  hand 
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high-speed  increases  the  consequence  of  the  ever-present  derailment  threat, 
road  maintenance  is  particularly  costly  if  freight  must  also  move  on  the 
high  quality  rails.  Magnetically  levitated  vehicles  of  either  the  attractive 
or  repulsive  type  are  less  susceptible  to  derailment  and  use  non-contact 
suspension,  but  seme  are  statically  unstable  and  all  face  the  cost  and  weight 
of  cryogenic  coil  cooling.  In  addition  unavoidable  electrical  discontinuity 
in  guideway  Joints  causes  field  interruptions. 

Air  cushion  vehicles  do  not  pose  such  problems.  They  are  statically 
stable  and  may  use  ordinary,  locally  rough,  non-conducting  guideway  materials 
such  as  concrete.  Derailment  is  not  a problem  a tracked  air  cusion  vehicle  is 
likely  to  face.  Low  loading  on  the  guideway  decreases  maintenance  costs.  Static 
cushion  vehicles,  however,  suffer  from  high  captation  drag,  induced  by  turning 
incoming  air  downward.  Even  this  problem,  though,  is  absent  in  a Tracked  Ram 
Air  Cushion  Vehicle  (TRACV).  Such  a vehicle  uses  dynamic  air  pressure  to 
generate  lift,  and  virtually  flies  in  ground  effect  down  its  guideway.  Its 
non-contact  suspension  is  not  as  noise-transmitting  as  a train,  and  the 
presence  of  the  ground  gives  it  far  better  lift/drag  characteristics  than 
aircraft . 

This  study  investigates  t'  ie  quality  of  a particular  TRACV  config- 
uration. As  a passenger  vehit  ossibility,  the  TRACV' s ride  quality  as 

measured  by  accelerations  and  suspension  travel  is  important.  A linear 
pitch-heave  model  of  the  vehicle  has  been  developed  and  the  longitudinal 
dynamics  are  analyzed  to  determine  the  vertical  accelerations  felt  by  the 
passengers  and  likelihood  of  vehicle -guideway  contact.  Three  configurations 
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are  considered:  l)  an  all-rigid  vehicle,  2)  a vehicle  with  moveable  control 
surfaces  suspended  by  springs  and  dampers,  and  3)  a vehicle  with  moveable 
control  surfaces  controlled  by  actuators  subject  to  active  feedback.  The 
basic  tradeoff  for  all  three  types  is  between  minimizing  vertical  accel- 
erations and  maintaining  adequate  clearance  between  the  vehicle  and  the 
guideway  without  unwieldy  control  power  requirements. 

PREVIOUS  WORK 

Since  the  TRACV  is  a fairly  recent  concept,  literature  discussing  its 
dynamics  is  not  abundant.  Dynamics  analyses  have  been  performed,  but  a 
• pitch- heave  analysis  which  accounts  for  the  considerable  coupling  has  not 
been  done.  Some  of  the  work  on  static  air  cushion  suspensions  is  indirectly 
applicable,  but  the  dynamics  are  considerably  different  and  results  do  not 
necessarily  correspond. 

T.  M.  Barrows  and  S.  E.  Widnall  (1970)  presented  a paper  which  laid 
out  the  aerodynamics  for  the  "ram  wing".  This  paper  discusses  ride  quality 
of  the  heave-only  constrained  vehicle  after  developing  expressions  for 
heave  and  pitch  forces  due  to  ground  disturbances. 

W.  E.  Fraize  and  T.  M.  Barrows  (1973)  performed  a system  definition 
study  which  hypothesized  a 154-foot  vehicle  with  a one  hundred  passenger 
capacity.  This  extensive  study  defined  all  aspects  of  a full-scale  system 
and  studied  feasibility  and  basic  dynamics.  Appendix  C of  this  reference 
discusses  a TRACV  suspended  on  one-degree-of-freedom  winglets  for  heave 
control  and  performs  a body  heave -win glet  dynamics  analysis  which  was  later 
found  to  contain  seme  errors.  Earlier  that  same  year,  Barrows  presented  a paper 
which  shows  that  almost  any  body  natural  frequency  and  damping  ratio  can  be 
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achieved  In  a suspended  TRACV  and  discusses  seme  active  control  concepts 
in  TRACV  suspension. 

P.  V.  Aidala  (197*0  studied  the  lateral  stability  of  the  TRACV  and 
determined  stability  derivatives.  His  report  also  compares  longitudinal 
forces  with  the  numerical  predictions  of  a one-dimensional  flow  mass 
conservation  theory  given  by  Boccardoro,  and  shows  that  lateral  and  longi- 
tudinal forces  may  be  considered  independent. 

H.  C.  Curtiss,  Jr.,  and  W.  F.  Putman  (1977)  present  a simplified  theory 
for  lift  and  moment  characteristics  and  show  good  agreement  between  this 
theory  and  wind  tunnel  and  towed  model  tests.  Their  report  pi-esents 
equations  for  lift  and  moment  stability  derivatives  and  describes  pitch 
and  heave  static  stability,  but  there  is  no  discussion  of  coupled  pitch- 
heave  dynamics. 

THE  VEHICLE 

The  vehicle  model  used  in  this  study  is  a full-scale  TRACV  with  a linear 
sloping  bottom  surface  and  haveing  the  cross-section  of  the  upper  half  of 
an  NACA  0021  airfoil.  It  is  150  feet  long  and  15  feet  wide  with  2.5-foot 
winglets.  The  guideway  is  12  feet  wide  and  about  7 feet  deep,  with  winglet 
running  lips  placed  at  a 45  degree  angle  to  horizontal.  The  size  of  the 
vehicle  (but  not  the  shape)  is  similar  to  that  of  the  Mitre  Corp.  Study 
(Fraize  and  Barrows,  1973).  Figure  1 shows  the  front  and  side  views  of 
the  vehicle.  The  bottom  slope  is  .028  rad  (about  1.5°)  and  the  bottom 
surface  clearance  at  the  trailing  edge  is  three  feet.  These  values  were 
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derived  from  5 = 1.4,  which  experiment  showed  to' give  good  L/D  character- 
istics and  good  correlation  between  theory  and  experiment.  The  venicle  is 
constrained  to  a constant  forward  velocity  of  300  fps  or  205  mph. 

The  winglets  ride  in  very  close  proximity  to  the  guideway  lips,  having 
a nominal  1.07  inch  clearance . With  the  vehicle  length  of  150  feet,  this 
close  clearance  will  necessitate  very  good  control  of  pitch  displacements; 
and  gap  variations  at  the  forward  and  trailing  corners  of  the  winglets  will 
prove  to  be  the  limiting  values.  The  reason  for  such  close  clearance  is 
that  the  winglets  essentially  act  as  a pressure  seal  to  keep  up  the  underside 
pressure.  Because  of  this  fact,  the  lift  variation  with  height  change  is 
strongly  dependent  upon  gap  variation;  the  stability  derivatives  support 
this  conclusion. 

OVERVIEW  OF  ANALYSIS 

This  thesis  first  reviews  the  various  terms  and  ratios  used  to  develop 
the  stability  derivatives.  The  latter  part  of  Chapter  H derives  the  control 
surface  derivatives  and  presents  the  basic  model.  It  then  describes  the 
guideway  roughness  modelling  process  and  describes  the  effect  of  size  of  the 
vehicle  on  short -wavelength  guideway  disturbances. 

Chapter  III  discusses  the  dynamics  of  the  rigid  vehicle  with  its  im- 
moveable winglets.  The  acceleration  spectral  density  is  compared  to  the 
Urban  Tracked  Air  Cushion  Vehicle  (UTACV)  standard,  and  ras  values  of  accel- 
eration and  winglet  gap  variation  are  presented.  These  cure  the  measures  by 
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which  this  study  determines  ride  quality  and  Chapters  III,  IV,  and  V follow 
similar  foraats. 

Chapter  IV  develops  the  passively  suspended  vehicle  using  massless 
wlnglets.  Linear  torsional  springs  and  dampers  are  attached  to  the  wing- 
lets  and  the  values  for  these  are  optimized  according  to  a criterion  based 
on  acceleration,  gap  variation,  and  winglet  deflection.  This  optimization 
procedure  is  a four-variable  gradient  search  performed  by  a digital  computer. 

Chapter  V presents  the  model  for  the  active  suspension.  Linear  optimal 
control  is  applied  to  find  the  optimal  feedback  gains  according  to  a per- 
formance index  similar  to  that  used  in  Chapter  IV.  Spectral  densities  are 
shown,  as  in  all  cases,  for  the  foremost  and  rearmost  passengers'  seats. 

The  feedback  gain  matrix  is  then  altered  in  a simple  check  to  find  the  effects 
of  reduced  sensor  capability  on  ride  quality. 
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CHAPTER  II 


MODELLING  THE  VEHICLE 

Given  that  the  TRACV  is  supported  by  aerodynamic  forces,  the  mathe- 
matical model  of  it  can  be  placed  in  a form  similar  in  some  respects  to 
those  of  aircraft.  Hence,  the  equations  of  motion  contain  coefficients 
of  lift  and  moment  as  well  as  stability  derivatives  similar  to  those  of 
aircraft.  At  this  point,  however,  most  of  the  similarities  end.  The 
parameters  of  interest  from  which  these  derivatives  are  calculated  are 
considerably  different  from  other  vehicles.  These  parameters  arise 
from  the  fact  that  the  TRACV  operates  in  ground  effect  and  they  are  ones 
which  can  be  measured  in  tests  and  expressed  explicitly  in  equations  of 
motion  and  derivatives.  They  are  also  useful  in  that  they  appear  in  non- 
dimensional  form  and  may  readily  be  applied  to  either  scale  models  or 
full-sized  vehicles.  In  addition,  the  axes  are  not  body  fixed,  but 
are  fixed  in  an  inertial  reference  and  translate  with  the  vehicle,  giving 
such  that  the  axes  are  not  free  to  rotate.  In  this  coordinate 
system,  the  word  "position”  refers  to  the  state  of  the  vehicle  in  pitch 
and  height  at  a given  point  in  time. 

This  Chapter  first  discusses  the  motion  variables  which  will  later 
be  used  as  state  variables,  then  explains  the  non-dimensional  ratio 
parameters  as  they  appear  in  the  work  of  Curtiss  and  Putman  (1977). 
Finally,  all  these  are  drawn  together  into  the  various  parts  of  the  model 
on  which  the  dynamics  analysis  is  to  be  done.  This  model  is  developed 
part  by  part  to  Include  attitude  stability  derivatives,  control  surface 


motions  and  guideway  roughness.  It  Is  hoped  that  this  development  will 
give  the  reader  a thorough  understanding  of  all  the  aspects  of  this 
unique  vehicle  being  studied. 

ATTITUDE  AND  RATE 

The  two  most  obvious  variables  to  present  themselves  are  the  height 
of  the  vehicle  and  its  pitch  attitude,  called  h and  0 respectively.  The 
height  variable,  h,  is  the  measure  of  the  vehicle's  height  in  feet  above 
some  reference,  the  value  of  which  is  divided  by  the  length  of  the  vehicle 
(called  chord  length,  or  c).  Thus,  h has  no  dimension.  Theta  (0)  is  the 
* pitch  attitude  of  the  vehicle  with  respect  to  some  reference,  measured  in 
radians . 

The  vertical  velocity,  6,  is  the  vertical  velocity  in  feet  per  second 
non-dimensionalized  by  the  forward  speed,  U.  This  variable  is  not  the 
same  as  the  time  derivative  of  h,  but  is  related  to  it  by  the  constant 
(£).  The  pitching  velocity,  q,  is  similarly  related  to  the  time  derivative 
of  9 by  the  same  constant. 

Thus,  there  are  four  variables  which  describe  the  state  of  the 
TRACV:  h,  0 , h,  q.  Upward  translation  and  nose  up  rotation  are  considered 
to  be  positive.  It  is  also  important  to  keep  track  of  the  relations 

£-§[&(»(]  »-§t&(e)l 

since  we  will  eventually  wish  to  use  common  variables  in  expressing  the 
equations  of  motion.  The  stability  derivatives  will  be  the  derivatives 
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Figure  1.  Geometry  of  or 


of  the  lift  and  moment  coefficients  with  respect  to  each  of  these  four. 


The  reference  against  which  each  of  these  motion  variables  is  to  be 
measured,  as  well  as  the  specific  motion  being  measured,  is  denoted  by- 
special  subscripting  or  prescripting.  For  instance,  the  height  of  the 
vehicle  center  of  gravity  above  some  inertial  reference  is  h^,  whereas 
the  relative  height  of  the  vehicle  c.g.  above  the  point  directly  below 
on  the  guideway  surface,  or  relative  height,  is  hr<  A full  description 
Of  various  notations  may  be  found  in  the  Nomenclature. 

RATIO  PARAMETERS 

In  keeping  with  the  idea  of  using  variables  which  may  be  conven- 
iently expressed  and  solved  for,  the  development  of  the  equations  of 
lift  and  moment  coefficient  are  functions  of  three  parameters:  a,  r , rx. 
These  three  are  ratios  of  different  areas  and  arise  from  the  fact  that 
the  lift  and  moment  equations  are  based  on  the  continuity  equation  for 
the  control  volume  which  is  the  airspace  between  the  vehicle  and  the 
guideway.  Air  enters  the  frontal  area  and  escapes  beneath  the  winglets 
and  out  the  trailing  edge  area  at  free  stream  pressure.  A detailed  des- 
cription of  the  theory  is  discussed  in  Curtiss  and  Putman  (1977 ).  Thus, 
the  most  convenient  way  to  express  lift  and  moment  coefficients  are  as  a 
function  of  the  ratios  of  entrance  and  exit  areas. 

The  first  ratio  is  <5,  which  is  defined  as  the  ratio  of  the  frontal 

area  less  the  exit  area  to  the  exit  area.  This  ratio  is  depicted  in 
Ai 

Figure  1 as  — . The  defining  relationship  is 


where  Ac  is  the  trailing  edge  exit  area,  W is  vehicle  width,  c is  chord 
length,  and  a is  the  slope  of  the  bottom  surface  with  respect  to  the 
guideway  expressed  in  radians . Note  that  for  small  angles  sin  a = a, 
so  hj  = coi.  In  Figure  1,  the  width  is  the  dimension  out  of  the  page. 

The  bottom  slopes  upward  toward  the  front,  giving  the  vehicle  positive 
angle- of -attack . 

The  second  parameter,  rQ,  is  a ratio  of  winglet  gap  area  to  area 

At  as  defined  in  Figure  1 when  the  winglets  are  parallel  to  the  guideway 

lip.  Its  defining  relation  is 

2 6c 
— o 

ro  - W a c 

The  term  6q  is  the  size  of  the  gap  between  the  winglet  tip  and  guideway 

lip  (see  Figure  2a)  and  the  factor  of  two  accounts  for  both  winglets. 

Hence,  r includes  the  total  winglet  gap  area,  not  just  that  for  one  side 
o 

Finally,  the  third  parameter  is  rx,  which  expresses  the  change  in 
winglet  gap  area  as  the  winglets  take  on  nonzero  angles  with  respect  to 
the  guideway  lip.  It  is  given  by 

-2  c 0 c 
ri  W a c 

and  represents  the  difference  in  winglet  gap  area  between  a non-parallel 
winglet  and  a parallel  one  divided  through  by  area  Ax  of  Figure  1.  Theta 
is  the  angle  the  winglets  make  with  respect  to  the  guideway,  or  simply 
the  body  pitch  angle  perturbation. 


Like  other  aerodynamic  vehicles,  the  TRACV  may  require  surfaces  of 

sane  sort  which  exert  controllable  forces  on  the  vehicle.  Since  the  motion 

variables  of  interest  are  pitch  and  heave  and  their  rates  of  change  with 

time,  surfaces  which  control  each  of  these  motion  separately  are  most  desirable 

fran  a control  standpoint.  Although  several  types  are  possible,  the  form 

which  appears  to  be  feasible,  realizeable,  and  giving  uncoupled  h and  ff  control 

is  a two-degree- of- freedom  winglet  rotation.  The  first  motion  is  rotation 

about  an  axis  parallel  to  the  vehicle  longitudinal  axis  and  is  designated 

6 ; the  second  is  a rotation  about  an  axis  perpendicular  to  the  longitu- 

wi 

dinal  axis  of  the  vehicle  located  at  about  a forty  percent  chord  position, 

designated  6 . These  motions  are  depicted  in  Figures  2a  and  2b  respectively. 

wa 

STABILITY  DERIVATIVES 

In  order  to  arrive  at  a reasonably  complete  and  yet  workable  model, 
we  must  make  a few  assumptions.  The  first  of  these  is  that  we  are  in- 
terested in  small  perturbations  from  the  steady-state  only.  Since  large 
motions  while  cruising  necessarily  mean  poor  ride . quality , they  must  be 
prevented  in  the  first  place  by  a tight  feedback  control  or  a suspension. 

The  second,  a direct  result  of  the  first,  is  that  the  lift  and  moment 
may  be  described  by  a Taylor  series  truncated  after  the  first  order  terms, 
since  second  order  and  higher  terms  are  vanishingly  small  in  a small 
perturbation  analysis. 

As  can  be  seen  in  Appendix  A,  the  steady-state  and  perturbation  parts 
of  the  problem  may  be  separated.  The  equations  of  motion  are  then  expressed 
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as  the  products  of  derivatives  and  corresponding  attitude  perturbations. 
These  derivatives  are  expressed  as  a series  of  partial  derivatives  with 
respect  to  each  of  the  three  area  ratio  parameters.  For  instance,  the 
change  in  lift  with  variation  in  height  is  expressed  as 


+ + 

ar.  • dr,  • ay- 


ah 3 a 3h  Uio  3h  1 3h 
and  in  similar  fashion  for  the  remaining  three  motion  quantities,  0,  h, 
and  q.  Each  of  these  partials  may  be  found  using  equations  found  in 
Appendix  B. 

In  addition  to  the  derivatives  with  respect  to  each  of  the  motion 
variables,  we  must  consider  the  effects  of  change  in  control  surface 
position  on  lift  and  moment.  There  are  four  such  terms:  lift  and 
moment  derivatives  with  respect  to  each  of  the  two  winglet  motions, 


and  6 


As  before,  these  are  expressed  as  a series  of  partial 


derivatives  with  respect  to  the  three  ratio  parameters;  for  example: 
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A full  development  of  these  is  also  fohnd  in  Appendix  B.  We  now  state 


the  third  assumption:  the  effect  of  winglet  motion  rate  (6^.  and  6’  ) 

w W ( 


on  vehicle  dynamics  is  negligible.  This  assumption  simplifies  the  model 
# 

and  seems  reasonable  since  the  damping  effects  of  winglet  rate  is 
probably  small  with  respect  to  steady-state  fqrees. 

The  result  is  a coupled  pair  of  fourth  order  equations  in  the  variables 


h,  q,  h,  and  0 with  the  control  variables  6 and  fi 


The  accelerations 
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of  the  vehicle  are  strictly  dependent  upon  the  relative  attitude  and 
velocity  of  the  vehicle  with  respect  to  the  guideway.  The  perturbation 
equations  of  motion  appear  in  the  following  form: 
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These  equations  are  further  refined  in  Appendix  A. 

MODELLING  TST1  GUIDEWAY 

Considering  that  this  is  a ride  quality  study,  the  effect  of  the 
guideway  roughness  on  vehicle  motions  is  a necessary  part  of  the  overall 
model.  Any  control  design  or  suspension  must  account  for  the  spectrum 
of  the  disturbances;  the  values  for  acceleration  and  suspension  displace- 
ment will  be  largely  dependent  upon  the  type  of  surface  over  which  the 
vehicle  will  travel.  As  the  entire  model  so  far  is  linear,  some  linear 
model  of  guideway  roughness  is  desirable. 

A search  of  the  literature  shows  that  a commonly  used  model  for 
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tt  x 10  ft.  is  a typical  value;  this  seems  a likely  value  for  a guideway 


constructed  for  the  TRACV  and  not  one  which,  as  a design  specification, 
would  require  strict  construction  tolerances,  driving  the  construction 


cost  out  of  practical  limits.  The  above  relationship  can  also  be  ex- 
pressed as  i = |h  (s)|  a.  It  then  follows  that  if  we  let  h (s)  = 

6 6 6 s 

(ignoring  phase  relationships)  then,  once  again  involving  the  small 


angle  assumption  (sin  6=6)  then  we  can  say 
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Thus,  it  seems  that  the  guideway  appears  as  white  noise  filtered  by  an 
integrator  in  height  and  simply  held  at  a constant  value  over  sill  frequencies 
in  pitch.  Using  white  noise  as  an  input  for  the  overall  model  simplifies 
the  analysis  and  control  problem  since  it  has  an  even  spectrum  and  does 
not  affect  closed  loop  characteristic  modes  of  motion.  Standard  methods 
for  control  such  as  found  in  Bryson  and  Ho  (1969)  may  be  directly  applied. 
FINITE  PAD  LENGTH  EFFECT 

The  above  guideway  model  is  incomplete,  however.  It  has  a drawback 
in  that  such  a definition  seems  to  say  that  at  any  instant  of  time  and 

0 

regardless  of  frequency,  the  entire  segment  of  guideway  under  the  vehicle 
has  a uniform  height  and  pitch.  For  instance,  a stairstep  in  the  guideway 
would  not  appear  to  the  vehicle  as  an  obstacle  to  climb,  but  as  a step 
change  (in  the  guideway  model)  in  the  height  of  the  entire  guideway.  If 
the  vehicle  bad  a point  contact  suspension  this  would  work,  but  in  a 
model  which  is  of  finite  length  and  attitude  and  rate  dependent,  the 
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result  is  infinite  rtns  accelerations  in  the  vehicle.  This  cannot  be,  but 
there  is  a correction  for  this  error  kncwn  as  the  finite  pad  length  effect. 

If  we  assume  that  that  which  primarily  affects  the  vehicle  is  the 
average  of  each  of  the  guideway  variables  over  the  length  of  the  vehicle, 
then  the  very  short  wavelength  disturbances  seem  to  be  filtered  out  by  the 
size  of  the  vehicle,  as  one  might  reasonably  expect.  This  assumption  is 
not  an  attempt  to  model  the  unsteady  aerodynamics  of  the  TRACV  which  are 
assumed  to  be  unimportant , and  is  made  by  way  of  analogy  to  other  levitated 
vehicles.  The  result  given  is  thoroughly  developed  in  Appendix  C (Riblich, 
Captain,  and  Richardson;  1967). 

The  effect  of  taking  the  average  of  guideway  variables  is  to  seemingly 

tsin 

— - — 1 . This  attenuation 

can  be  closely  approximated  by  the  magnitude  function  of  a second  order 
filter,  such  as  the  one  depicted  in  Figure  3.  A second  order  filter  was 
chosen  because  it  had  sufficient  dropoff  to  insure  that  there  were  no 
infinite  acceleration  in  the  total  model  (as  would  have  been  the  case  with 
a first  order  model)  although  it  slightly  overestimates  the  attenuation 
of  the  I f>unction*  This  filter  has  the  transfer  function  (in  real 
time) 


The  linear  model  of  the  guideway  as  it  appears  to  the  vehicle,  then, 
la  white  noise  filtered  by  an  integrator  and  further  filtered  by  a second 
order  finite  pad  length  approximation  for  height  disturbances,  and  the 
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time  derivative  of  the  same  for  pitch  disturbances.  The  relationship  is 
illustrated  in  block  diagram  form  in  Figure  4. 


Figure  4.  Block  Diagram  Representation  of  Guideway 
Model. 


The  advantage  of  using  such  a linear  model  is  that  it  does  not  unduly 
complicate  the  understanding  of  vehicle  dynamics  and  it  does  not  affect 
a Gaussian  process  (white  noise)  other  than  in  magnitude.  The  guideway 
spectra  may  then  be  expressed  by  the  following  differential  equations. 

9 _ + 13.7  8 + 45.3  9 Tl(t) 

6 o o ^ 


where  T|(t)  is  a white  noise  source  with  unity  spectral  density.  This  is 
a third  order  set  of  equations  which  completely  describes  the  guideway 
variations  as  seen  by  the  vehicle  traversing  over  it. 


TOTAL  MODEL 


The  total  model  is  a fourther  order  set  of  vehicle  equations  controlled 
by  two  inputs  and  disturbed  by  a third  order  white  noise  filter  which  models 
the  guideway.  It  is  expressed  in  state  variable  form  as  x = Ax  + Bu  + Cw 
where  A,  B,  and  C are  matrices  and  x,  u,  w are  vectors.  These  are  presented 
in  Appendix  D.  The  stability  derivatives  are  derived  from  lift  and  moment 
coefficients  and  a first  order  Taylor  series  assumption.  Since  we  are 
interested  only  in  small  perturbations  about  some  steady-state,  we  can  elimin 
ate  the  steady-state  terms  from  the  equations  of  motion.  The  guideway  is 
assumed  to  have  the  spectrum  of  a typical  smooth  highway  and  the  effect 
of  vehicle  size  on  small  bumps  has  been  approximated  by  assuming  that  the 
average  value  of  the  guideway  under  the  vehicle  which  is  of  interst. 

The  complete  model  as  presented  does  have  its  drawbacks.  The 
linearity  of  the  system  will  not  predict  the  changes  in  lift  coefficient 
derivatives,  for  instance,  when  the  gap  is  altered  drastically  in  violation 
of  the  small  perturbation  assumption.  However,  the  theory  shows  that  the 
derivatives  remain  approximately  constant  over  gap  height  change.  For 
example,  the  steady-state  gap  of  one  inch  presently  used  corresponds  to 
rQ  = 0.5.  Varying  rQ  as  much  as  ± 0.475  changed  none  of  the  lift  deriv- 
atives by  more  than  three  persent.  All  the  derivatives  were  within  a 

&CM 

fifteen  percent  tolerance  except  for  which  varied  by  thirty-five 

dh 

percent.  This  seems  to  say  that  a linear  assumption  is  reasonably  good. 

Also,  there  is  no  provision  for  air  turbulence;  the  model  assumes  a 


21 


stationary  air  mass  in  the  guideway.  Though  this  last  assumption  is  not 
likely  to  be  the  case,  it  does  provide  a good  starting  point  for  the 
TRACV  dynamic  analysis  and  will  yield  a good  impression  of  what  the  dominant 
modes  of  motion  are  likely  to  be. 

I 4 
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Having  assembled  a complete  model,  the  first  item  of  business  is  to 
explore  the  dynamics  of  the  rigid  vehicle,  one  in  which  the  winglets  are 

, 

firmly  fixed  with  respect  to  the  body.  Doing  so  provides  an  initial  idea 


of  the  fundamental  vehicle  dynamics  and  whether  some  form  of  suspension 
is  necessary.  Setting  up  the  tools  for  this  analysis  also  paves  the  way 
for  later  more  complicated  analysis. 

EIGENVALUES  AND  EIGENVECTORS 

We  start  with  a look  at  the  eigenvalues  and  eigenvectors,  for  they 
describe  the  vehicle  characteristics  in  detail.  Table  I gives 
these  for  the  rigid  vehicle.  From  these  values  we  can  see  there  are  two 
oscillatory  modes  in  the  vehicle:  one  low  frequency,  very  lightly  damped 
mode  and  a second  higher  frequency,  moderately  damped  mode.  Because  there 
are  two  complex  pairs,  only  two  eigenvectors  need  be  expressed;  the  other 
two  are  the  complex  conjugates,  of  these.  The  pitch  attitude  vari- 
able (0)  is  the  base  for  each  eigenvector. 

The  rigid  vehicle  has  two  natural  frequencies  which  are  fairly  close 
together,  neither  of  which  are  strongly  damped.  The  eigenvectors  indicate 
that  the  height  variation  and  pitch  variation  are  nearly  in  phase  for  the 
low  frequency  mode.  This  means  that  the  nose  comes  up  with  the  rising  of 
the  c.g.  in  the  guideway  such  that  the  front  end  of  the  vehicle  has  the 
largest  excursions.  One  might  infer  that  the  winglet  gap  at  the  front 
will  have  the  greatest  variation,  since  forward  of  the  c.g.,  the  pitch 
and  heave  oscillations  add  to  each  other.  The  attitude  mode  ratio  for 
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Figure  1.  Transient  Response  of  Rigid  Vehicle. 


the  low  frequency  eigenvector  is  ^ = 0. ^b,  indicating  that  a 0.1  radian 
pitch  variation  coincides  with  an  8.1  foot  variation  in  c.g.  height.  For 
the  150  foot-long  vehicle,  this  motion  appears  as  a rotation  about  an 
effective  center  at  about  the  percent  chord  position,  almost  at  the 
tail  end.  Figure  1,  a time-domain  step  response,  shows  this  in-phase 
relationship.  Having  an  effective  center  of  rotation  at  .9^c  suggests 
that  the  forward-most  seat  in  the  cabin  will  experience  the  worst  acceler- 
ations, having  the  longest  moment  arm. 

On  the  other  hand,  the  high  frequency  mode  is  sixty  degrees  out  of 
phase  and  has  a mode  ratio  g-  which  is  six  times  greater  than  that  of  the 
low  frequency.  This  will  mean  that  this  frequency  appears  as  nearly  all 
heave  and  the  small  pitch  component  does  not  add  as  much  vertical  acceler- 
ation to  any  vehicle  location  as  does  the  low  frequency  mode.  This  effect 
is  apparent  in  the  spectral  density  plots  which  follow,  particularly  if 
one  remembers  that  part  of  the  high  frequency  reduction  at  .90c  is  due 
to  the  reduced  reinforcement  of  the  nearby  low  frequency  mode. 

SPECTRAL  DENSITY  AND  RIDE  QUALITY 

Table  II  shows  the  A and  C matrices  for  the  rigid  vehicle.  Since 
there  are  no  control  motions  (winglet  movements),  the  B matrix  is  of  no 
concern.  Using  these  matrices  and  Cramer's  rule,  the  transfer  functions 
may  be  found  for  the  various  state  variables  with  respect  to  the  noise 
vector  which  has  been  normalized  to  one.  The  magnitude  squared  of  this 
function  plotted  against  frequency  is  the  spectral  density.  For  vehicle 
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TABLE  2- 

OEHERALIZED  MATRIX  FORM  OF  RIGID  VEHICLE  STATE  EQUATIONS 


accelerations,  we  use  |s  h^|a  as  our  c.g.  spectral  density,  and  a linear 
combination  of  hj  and  0^  will  give  the  acceleration  transfer  function 
at  other  vehicle  locations.  These  spectral  densities  may  be  plotted 
against  a ride  quality  standard  which  will  show  at  what  frequencies  the 
accelerations  are  unacceptable. 

One  such  standard  is  the  Urban  Tracked  Air  Cushion  Vehicle  (UTACV) 
standard  (DOT  Specification,  1972)  which  sets  an  upper  limit  on  the 

•» 

magnitude  of  the  spectral  density  plot.  Figure  2 shows  the  spectral 
density  of  accelerations  at  the  c.g.  as  compared  to  the  UTACV  limit. 

It  shows  that  the  accelerations  in  the  range  from  0.764  - 2.24  Hz  are 
unacceptable.  At  the  rearmost  seat,  near  the  apparent  low  frequency 
mode  center  of  rotation,  we  would  expect  that  the  spectral  density  plot 
would  show  a diminished  low  frequency  peak  and  overall  reduced  accelera- 
tion according  to  the  prediction  of  the  eigenvectors.  Figure  3 shows  this 
to  be  the  case,  although  the  standard  is  scarcely  met. . These  figures 
and  similar  ones  in  later  paragraphs  are  in  non-dimensional  frequency. 

The  defining  relationship  is  1 Hz  - 9*15  rad/r.  The  vertical  axis  has 
units  of  db  (20  log  ). 

RMS  ACCELERATION  AND  GAP  VARIATION 

In  addition  to  spectral  density  comparisons,  an  important  indicator 
of  ride  quality  is  the  root  mean  squared  (RMS)  acceleration.  Whereas  the 
spectral  density  plot  shows  acceleration  level  at  a frequency,  the  rms 
value  gives  the  overall  acceleration  felt  by  the  passenger.  Both  indicators 
are  necessary  for  a complete  picture.  If  G(s)  is  the  transfer  function  of 
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interest,  then  the  mean  squared  value  is  given  by 


E(xa)  = J*  G(s)  G(-s)  ds;  s = j (0 

This  integral  has  been  solved  in  closed  form  and  tabulated  for  systems  up 
to  tenth  order  (Newton,  Gould,  and  Kaiser;  1957)»  At  the  c.g.,  the  rms 
acceleration  was  found  to  be  0.121  g,  a fairly  rough  ride.  The  rearmost 
seat,  at  .90c,  proved  to  be  far  better  as  predicted,  having  rms  accelera- 
tion of  0.0055  g. 

Not  only  are  we  interested  in  the  rms  value  of  the  accelerations,  but 
the  rms  variation  in  the  gap  at  each  end  will  give  an  idea  of  frequency  of 
contact  with  the  guideway.  Once  again,  the  Gaussian  white  noise  assumption 
is  useful,  because  the  gap  variation  is  Gaussian  as  well,  a direct  result 
of  having  a linear  gap  to  white  noise  transfer  function.  The  rms  value  is 
related  to  the  standard  deviation  of  the  Gaussian  distribution,  so  a 
knowledge  of  the  steady-state  gap  size  and  the  rms  variation  will  indicate 
what  percentage  of  the  time  the  winglets  are  in  contact  with  the  guideway. 

The  rigid  vehicle  leading  edge  rms  gap  variation  is  14.7  inches  and 
the  trailing  edge  rms  gap  variation  is  6.07  inches  as  compared  to  a 
steady-state  value  of  about  one  inch.  Since  the  vehicle  can  not  appear 
below  the  guideway,  these  values  tell  us  that  the  vehicle  is  likely  to 
spend  much  of  its  time  in  contact  with  the  guideway,  or  simply  put,  it 
will  "bottom  out"  frequently,  particularly  at  the  front,  primarily  be- 
cause of  the  influence  of  the  very  lightly  damped  low  frequency  mode. 

The  reader  should  remember  that  these  numbers  are  the  result  of  an  un- 
constrained linear  model  and  while  they  are  quite  good  as  estimates  of 
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what  will  occur,  they  cannot  exactly  describe  the  actual  nonlinear  problem 
of  guideway  contact. 


=1 


The  conclusion  to  be  drawn  from  this  rigid  case  analysis  is  that  from 
»n  measures  the  ride  quality  of  the  rigid  vehicle  is  unacceptable.  There 
are  frequencies  in  all  positions  of  seats  where  the  accelerations  are  un- 
comfortable. The-  overall  ride  is  rough,  and  the  vehicle  tends  frequently 
to  bottom  out.  The  vehicle  has  one  oscillatory  mode  which  is  very  lightly 
damped  at  about  1 Hz  which  is  in  large  measure  the  cause  of  these  problems, 
only  to  be  aided  by  another  mode  of  light  damping  nearby  at  1.5  Hz.  It  is 
readily  apparent  that  some  form  of  suspension  or  control  is  necessary  if 
the  TRACV  is  to  be  implemented.  The  Chapters  that  follow  explore  this 
possibility  in  both  the  passive  and  active  sense  and  compare  those  cases 
to  this  one  to  determine  what  degree  of  improvement  can  be  made. 


Since  it  is  clear  that  the  rigid  vehicle  is  unacceptable,  the  ways 
in  which  the  ride  quality  can  be  made  acceptable  become  important.  If 
the  winglets  are  given  freedom  to  move  about,  then  they  can  be  used  to 
control  winglet  gap  and  thus  vehicle  motions,  since  the  theory  of  lift 
and  moment  indicates  that  the  forces  are  very  strongly  a function  of 
the  winglet  gap.  The  motions  of  the  winglets  as  described  in  Chapter  II 
can  be  controlled  by  two  means:  passive  suspension  and  active  control. 

In  this  chapter  the  former  case  is  considered;  a model  is  developed  and 
the  ride  quality  of  this  model  is  examined.  Since  the  passive  suspension 
is  much  less  expensive  to  implement  than  active  control,  if  their  per- 
formances can  be  made  nearly  equal,  then  passive  suspension  is  preferable. 
THE  REVISED  MODEL 

The  passively  suspended  model  can  be  considered  an  uncontrollable 
dynamic  system  as  opposed  to  a system  with  control  surfaces,  with  the 
winglets  restrained  by  springs  and  dampers.  We  assume  massless  winglets 
for  simplicity  because  they  will  be  very  small  in  mass  compared  to  the 
vehicle  and  their  natural  frequencies  will  be  so  high  that  for  the  present 
purpose  they,  are  infinitely  fast.  There  are  then  two  first  order  equations 
which  a*s  added  to  the  equations  of  the  rigid  case,  one  for  each  degree-of- 


freedcm  of  the  winglets.  The  motions  of  the  winglets  A6  and  A6  are 

wx  w, 

considered  state  variables  and  the  coefficients  become  part  of  the  A 
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matrix.  After  finding  values  of  spring  and  damper  constants  which  give 
the  best  performance,  the  analysis  of  such  performance  parallels  that 
of  Chapter  III. 

We  start  by  writing  a force  balance  equation  for  each  winglet 
motion.  Once  again,  only  perturbations  are  of  present  concern.  The 
hinge  moment  on  the  winglet  for  each  mode  of  motion  is  related  to  the 
pressure  distribution  under  the  vehicle  and  is  a function  of  the  motion 
and  control  variables. 

Bi  A6Wj  + Ka  A6Wi  = HM  (&,  q,  h,  0 , 6^,  6^) 


Ba  + Ka  A6w  = HM  (fi,  q,  h,  0,  6.6) 


In  the  above  equations,  HM  and  HM  refer  to  the  hinge  moment  around  the 
first  and  second  winglet  axes,  respectively.  For  the  first  of  these  two, 
the  hinge  moment  is  related  to  the  lift  since  the  pressure  generating 
lift  places  a moment  on  the  winglet. 
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We  note  that  x is  the  distance  from  the  axis  to  the  center  of  pressure 
cp 

cn  the  winglet  in  feet  and  Aw  is  the  total  area  of  the  winglets  (including 

the  area  of  both  winglets).  Since  the  theory  assumes  that  the  pressure  at 

the  gap  is  that  of  free  stream  and  the  pressure  near  the  root  of  the  winglet 

is  that  under  the  vehicle,  there  must  be  some  gradient  of  pressure  along  the 

width  of  the  winglet.  Assuming  a linear  gradient,  as  approximation  the 

W 

center  of  pressure  is  at  or  0.833  feet.  This  value  yields  P =*  33^+0  ft. -lbs 
The  second  type  of  hinge  moment  is  also  related  to  previous  lift  and 
moment  coefficients , but  in  a different  way. ' 

~ _ . \r  L Aw  M 

HM  = (*c<r#  “ *0)  c ~2  s + ~2  S 

In  the  above  equation,  the  first  term  describes  the  pressure  at  the  c.g. 
due  to  lift  acting  on  the  winglet,  the  second  relates  the  effect  of  the 
lengthwise  pressure  distribution  which  causes  the  aerodynamic  moment  on 
the  vehicle.  Again  calling  on  the  assumption  that  *c  r “ xQ>  the  first 
term  drops  out  and  we  get 
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Because  of  our  earlier  notation  CL  = — — the  radius  of  gyration 

term  must  be  carefully  reinserted  to  insure  consistency.  This  development 
yields  R » 3753  ft. -lbs. 

The  basic  equations  of  motion  in  state  variable  form  appear  as 
i = Ax  + Cw  where  A,  instead  of  being  a seventh  order  matrix,  is  a ninth 
order  matrix.  The  eighth  and  ninth  rows  come  from  the  equations  below 
which  are  once  again  written  in  non-dimensional  time. 

' -rrCT.  Ah'  + ^CT  A6'  + P CT  Ah  + P CL  A0 
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Notice  the  change  in  notation  to  the  new  damping  constants,  Bx  and  Ba. 
This  is  to  indicate  that  these  two  have  also  been  shifted  to  non-dimen- 
sional time.  The  full  state  matrix  form  with  the  two  added  equations 
be  found  in  Appendix  D. 
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OPTIMIZING  THE  PASSIVE  SUSPENSION 

The  ride  quality  of  the  passively  suspended  TRACV  is  dependent  on 
the  four  spring  and  damper  constants;  Kx,  Ka,  B1,  and  Bs.  The  values 
which  give  the  best  possible  ride  quality  must  be  determined,  but  the 
process  involved  is  not  very  simple  because  of  a ninth-order  model  and 
a ride  quality  which  is  a function  of  four  variables.  The  most  obvious 
method  of  solution  is  a gradient  search  on  a computer,  since  the  analyt- 
ical method  would  involve  multiple  determinants  of  a ninth  order  system 
and  as  a result  the  likelihood  of  multiple  algebraic  errors.  A gradient 
search  which  simply  looks  for  the  direction  of  largest  improvement  in  a 
cost  function  and  adjusts  the  four  variables  accordingly  would  be  fast 
and  relatively  simple  to  program. 

The  routine  which  was  used  was  the  simplest  gradient  method  possible, 
based  on  the  procedure 


k + 1 


- « ( 


J*-  J* 
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Where  x is  the  particular  variable  being  incremented,  k is  the  iteration 
n 


counter,  c in  the  step  size  constant  for  x , and  J is  the  cost  function. 
1 n n 

The  procedure  is  to  find  the  partial  derivative  of  the  cost  function  with 
respect  to  each  of  the  four  variables,  adjust  by  the  corresponding  «Q, 
and  use  the  result  as  the  stepsize  for  each  variable. 

The  cost  function  which  was  used  placed  a cost  on  six  values.  These 
values  are  the  rms  acceleration  in  heave  and  pitch,  rms  gap  variation 


found  for  a given  A matrix  by  using  Cramer's  rule  to  find  the  transfer 
function  of  each  with  respect  to  unity  white  noise  and  calculating  the 
mean  squared  error  integral  of  that  transfer  function.  The  value  of 

V _ 1 

was  taken  as  five  percent  larger  than  rather  than  the  actual 
previous  value  to  insure  that  a large  previous  step  does  not  adversely 
influence  the  derivatives  at  the  present  point.  Computationally,  using 
old  values  would  have  been  less  expensive,  but  the  cost  margin  was  deemed 
small.  The  cost  function  used  is  described  by 

1 T 

J = ^ I x Zx  dt 

2 o 

where  x is  a vector  of  the  state  variables  and  Z is  a diagonal  matrix  of 
weighting  constants.  The  values  of  Z were  set  by 

Z(»,n)  - 

max 

which  is  the  inverse  of  the  desired  maximum  value  of  each  of  the  respective 

state  variables.  The  values  chosen  were  h*  2 .Olg,  Ah-.  = 0.33  inches. 

max  Gap 

and  no  weighting  on  control  deflections  or  pitching  accelerations.  These 
values  were  actually  the  result  of  first  having  solved  the  active  optimal 
control  problem  and  finding  out  what  set  of  weighting  coefficients  worked 

best. 

Using  this  cost  function,  an  algorithm  for  the  optimization  of  the 
spring  and  damper  constants  was  developed.  Figure  1 gives  a flow  diagram 
for  the  algorithm.  When  implemented  on  the  IBM  370  computer  using  APL, 
the  method  converged  in  approximately  twenty  steps.  The  cost  function 

.1 
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had  severed,  local  mini mums  in  the  same  general  vincinity,  so  the  one  with 
the  smallest  rms  vertical  acceleration  was  chosen.  Only  one  cost  function 
was  used;  it  is  possible  to  find  different  optimal  values  by  adjusting  the 
weighting  factors. 

With  the  cost  function  used,  the  solution  of  the  optimization  problem 


is  as  follows: 


Kx  - 3472 


Bj  = 1652 


ft. -lbs. 


ft. -lbs. 
rad/sec 


Ka  = 374.1 


ft. -lbs. 


Ba  = 1115 


ft. -lbs. 
rad/ sec 


Since  the  axis  of  rotation  of  A6  is  at  the  c.g.,  the  lift  causes  no 

w8 

steady-state  moment  on  Ka  and  the  fact  that  there  is  no  steady-state  moment 

to  counter,  the  spring  Ka  has  a zero  steady-state.  Spring  Kx,  on  the  other 

band,  must  supply  a steady-state  moment  to  the  winglet  to  counter  the  effect 

of  lift.  This  value  is  given  by 

K,  = HM  = P.C_  = 15225  ft. -lbs. 

0 0 Lo 

Since  these  are  the  total  spring  and  damper  values ; they  may  spread  out 
over  several  springs  and  dampers  for  each  motion  if  necessary.  This  set 
of  results  seems  to  be  the  optimum  passive  suspension  for  the  given  set  of 
weighting  coefficients.  It  was  felt  that  a better  suspension  might  be  found 
by  repeating  the  entire  process  with  several  sets  of  weighting  parameters. 
Rather  than  grind  away  numerically  at  considerable  costs,  however,  a 
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local  search  around  this  optimum  was  made  by  using  educated  guessing,  and 
no  better  tradeoff  between  vertical  acceleration  and  gap  variation  could 
be  found. 

With  the  spring  and  damper  constants  determined,  the  overall  model 
is  easily  described  in  state  variable  form  with  two  more  variables  than 
in  the  rigid  case;  A6w  and  A6w  are  becoming  state  variables.  We  now 
have  a ninth  order  system  which  may  be  evaluated  using  the  same  criteria 
as  used  in  Chapter  III. 

EIGENVALUES  AND  EIGENVECTORS 

Once  again  we  start  with  a study  of  the  eigenvalues  and  eigenvectors. 
For  the  passive  suspension  they  differ  greatly  from  the  rigid  case,  as 
shown  in  Table  I.  The  eigenvalues  show  a marked  change  from  very  low 
damping  to  moderate  damping  and  drastic  lowering  of  natural  frequencies. 
For  the  rigid  vehicles,  the  two  frequencies  occurred  at  two  and  three 
hertz;  in  the  passive  case  they  occur  at  0.036  and  O.325  Hz,  so  that  they 
are  not  only  lower,  but  they  are  relatively  further  apart.  Spreading 
them  apart  and  increasing  the  damping  means  that,  on  the  whole,  the  accel- 
eration will  be  much  lower  and  more  nearly  uniform  over  the  range  of  fre- 
quencies . 

The  eigenvectors  also  show  a marked  change.  The  mode  ratios  of 
heave  to  pitch  for  the  two  frequencies  are  much  larger,  indicating  that 
the  degree  of  coupling  is  drastically  reduced.  The  high  frequency  mode 
has  a ^ mode  ratio  which  is  sixty-four  times  greater  than  the  ratio  for 
the  rigid  case,  implying  that  the  high  frequency  mode  is  nearly  all  heave. 
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TABLE  1 

EIGENVALUES  AND  EIGENVECTORS  OF  PASSIVE  SUSPENSION 


Eigenvalues 
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-5.5 

0 

-.967  + i 1.79 
-.967  - 1 1.79 
-.077  + i .215 
-.077  - 1 .215 

Winglet  Modes 
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Eigenvalues  in  units  of  sec 
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PASSIVE  CASE  RESPONSE  TO  HDOT*~  1 AND  6 


PLOT  OF  HEIGHT  AND  PITCHl  HEIGHT  IS  SMOOTH  CURVE 


Figure  2.  Transient  Response  of  Passively  Suspended  Vehicle 
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For  the  low  frequency  mode,  the  opposite  is  true,  the  — mode  ratio  for  the 

h 

passive  case  is  over  ten  times  greater  than  that  of  the  rigid  case  and  the 
low  frequency  mode  is  nearly  all  pitch  motion. 

These  values  suggest  that  the  motions  of  A6  act  so  as  to  decouple 

wi 

almost  entirely  the  two  modes  of  motion  of  the  body  while  they  move  to 
track  variations  in  the  guideway,  keeping  the  vehicle  body  as  motionless 
as  possible.  The  high  frequency  damping  ratio  has  increased  from  the  rigid 
case  value  of  0.23  to  0.48  and  the  low  frequency  ratio  has  improved  from  the 
rigid  case  0.08  to  0.34,  indicating  that  these  winglet  motions  also  reduce 
the  number  of  oscillations  in  the  vehicle  body.  Because  of  the  very  low 
frequency  of  the  pitch  mode,  the  accelerations  caused  by  the  pressure  of 
this  mode  are  apt  to  be  snail  (given  that  the  excursions  are  not  very  large) 
and  hence  the  accelerations  felt  by  the  passengers  will  be  mostly  due  to 
heave  and  should  be  nearly  equivalent  everywhere.  Figure  2 illustrates 
a typical  time-domain  step  response  of  the  passively  suspended  vehicle. 
ACCELERATION  SPECTRAL  DENSITY 

Shown  in  Figures  3 and  4 are  the  vertical  acceleration  spectral  den- 
sities at  the  c.g.  and  .90c  positions,  respectively,  which  show  a definite 
improvement  over  the  rigid  case.  Not  only  does  the  curve  meet  the  UTACV 

specification,  but  it  is  better  by  at  least  a factor  of  twenty  at  all 
# 

frequencies.  There  are  no  sharp  peaks,  indicating  that  the  damping  is 
Improved  as  well.  These  figures  would  suggest  a very  acceptable  ride 
quality  as  a result  of  the  passive  suspension  by  the  low  level  of  acceler- 
ation at  all  frequencies. 
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Passive  Suspension  Acceleration  Spectral  Density  at  .40c 
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RMS  ACCELERATION.  GAP  VARIATION.  AND  CONTROL  DEFLECTION 

The  acceleration  at  the  c.g.  was  greatly  improved  by  the  passive 
suspension.  In  the  rigid  vehicle  the  rms  vertical  acceleration  was  found 
to  be  0.121g,  while  the  passive  suspension  improved  this  to  0.0072g,  a 
very  smooth  ride.  The  rearmost  seat  was  about  the  same  as  the  c.g.  in 
the  passive  suspension,  as  predicted,  experiencing  0.009g  accelerations 
which  are  slightly  worse  than  the  rigid  case.  The  very  smooth  .90c 
acceleration  level  for  the  rigid  case  was  due  to  the  in-phase  coupling 
of  the  modes  which  is  destroyed  in  the  passive  suspension.  Hence,  the 
slight  loss  of  ride  quality  at  that  location  is  understandable. 

Another  significant  improvement  is  in  the  gap  variation.  The  leading 
edge  rms  gap  variation  is  one  percent  of  and  the  trailing  edge  value 
is  3.5$  of  the  rigid  case  value.  These  new  values  are  Ah^  * .143 
inches  and  Ah,^,  = .224  inches.  Since  the  steady-state  gap  is 

slightly  over  one  inch,  one  can  safely  say  that  the  winglets  rarely  touch 
the  guideway.  In  fact,  it  seems  that  there  is  seme  margin  in  these  values, 
and  that  adjusting  the  cost  function  to  weight  accelerations  more  heavily 
might  result  in  seme  set  of  springs  and  dampers  which  enlarge  this  variation 
slightly  and  yielded  still  lower  acceleration  levels.  It  is  interesting  that 
the  trailing  edge  has  the  larger  variation,  which  is  not  true  of  the  rigid 
vehicle.  It  is  clear  from  the  above  values,  however,  that  the  difference 
is  smaller,  as  a result  of  mode  decoupling.  The  decoupling  and  phasing  of 
h and  6 motions  is  further  illustrated  by  this  fact. 
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There  is  one  other  result  of  interest,  and  that  is  the  rms  deflection 


of  the  two  wlnglet  modes.  Thest  values  are  found  in  the  same  way  as  the 


others  and  give  an  idea  of  the  typical  size  of  the  control  surface  displace- 
ments. For  the  passive  suspension,  rms  A 6 is  0.0154  rad.  and  rms  A6 

wi  wa 

is  0.000534  rad.  The  typical  heave  control  is  small,  varying  through  a 


typical  ± 1°.  A deflection  of  one  degree  represents  a change  of  gap  on 
the  order  of  one -half  inch,  however,  and  that  is  significant.  This  seems 
to  tell  us  that  the  motions  of  the  winglet  are  such  that  the  heave  in  the 
guideway  is  countered  by  similar  motions  of  the  winglet  so  that  the  gap 


remains  constant  along  with  the  c.g.  while  the  winglet  and  guideway  move 


together.  The  typical  pitch  deflection  is  also  quite  small,  for  we  find 
that  the  value  of  0.000534  rad.  corresponds  to  .03°,  and  a gap  change  at 
the  vehicle  rear  of  about  0.6  inches.  Once  again,  the  winglet  deflection 
is  larger  than  the  gap  variation,  indicating  that  the  second  mode  of  motion 
is  also  probably  necessary  to  keep  the  vehicle  body  unaccelerated  and  to 
prevent  wingle t - guideway  contact. 

To  summarize,  we  have  built  a passive  suspension  by  allowing  the 
vinglets  two-degrees-of-freedom  but  restraining  each  mode  with  a spring 
and  a damper.  The  effect  of  this  configuration  (if  one  finds  optimal 


springs  and  dampers)  is  to  lower  the  frequencies  of  both  modes  while 


spreading  them  apart,  increasing  damping,  and  decoupling  the  modes  of  motion. 
In  terms  of  the  ride  quality,  the  new  configuration  is  very  good  and  a 
tremendous  improvement  over  the  rigid  vehicle.  The  overall  acceleration 
level  is  quite  low  and  there  are  no  particular  "bone- jarring"  frequencies. 
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: i 


: 


The  rate  of  vlnglet  contact  goes  from  almost  continuous  to  almost  never, 
and  the  size  of  springs  and  dampers  required  is  reasonable  and  permit 
only  very  small  control  deflections,  so  that  the  passive  suspension  does 
not  attempt  to  violate  the  assumptions  from  which  the  model  arose.  All 
in  all,  the  general  outlook  is  very  good  for  this  form  of  control. 
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CHAPTER  V 


ACTIVE  SUSPENSION 


The  active  suspension  is  one  in  which  the  winglets  are  controlled 
by  actuators  which  receive  signals  from  feedback  circuits  originating 
at  vehicle  position  and  rate  sensors.  In  this  chapter  the  linear 
optimal  control  problem  is  solved  using  a performance  index  very  similar 
to  the  one  in  the  previous  chapter.  The  form  of  the  equations  is  varied 
slightly  to  accommodate  the  solution  technique  in  that  the  winglet  motion 
equations  are  dropped  as  a result  of  assuming  that  the  winglet-servo 
response  to  control  inputs  is  of  sufficient  bandwidth  to  be  considered 
infinitely  fast.  The  method  of  solution  involves  the  solution  of  the 
matrix  Ricatti  equation  and  using  this  to  find  the  feedback  gains.  Solving 
this  problem  will  give  seme  idea  of  the  best  which  this  particular  vehicle 
can  possibly  do,  and  provides  an  index  against  which  other  types  of  sus- 
pensions may  be  compared.  Having  found  the  optimal  feedbacks,  we  will 
proceed  to  analyze  the  ride  quality  along  the  lines  of  previous  chapters. 
ACTIVE  SUSPENSION  MODEL 

This  model  assumes  a controllable  dynamic  system  with  white  noise 
disturbance.  The  state  equations  written  as  matrices  have  the  form 
£ = Ax  + Bu  + Cw  where  the  A and  C matrices  are  identical  to  those  of 
the  rigid  vehicle  model,  and  the  B matrix  is  the  set  of  derivatives  of 
lift  and  moment  coefficient  with  respect  to  each  of  the  two  control 
motions.  The  result  is  a seventh  order  system  with  two  control  inputs 
and.  one  disturbance  input. 


B i I 


Since  there  are  seven  state  variables,  this  formulation  assumes 
seven  sensors  from  which  the  feedback  signals  originate.  This  may  not 
always  be  feasible,  so  a comparison  is  done  later  with  a system  in  which 
some  of  the  sensors  are  turned  off.  The  seven  states  are  vehicle  heave 
rate  and  pitch  rate,  relative  heave  and  pitch,  and  guideway  heave  rate, 
pitch  rate,  and  pitch.  The  first  two  may  easily  be  sensed  with  on-board 
accelerometers,  the  output  of  which  is  integrated.  The  second  pair  could 
be  sensed  by  measuring  capacitance,  using  small  feelers,  or  even  by  care- 
fully measuring  under-vehicle  pressure  at  several  points  and  deriving 
attitude  from  the  lift  and  moment  relationships.  Of  these  possibilities, 
use  of  feelers  is  the  least  desirable  since  the  feelers  create  noise  and 
would  continually  require  replacement  because  of  wear. 

Infrared  ranging  of  the  relative  distance  is  difficult  because  the 
track  is  subject  to  differential  heating  which  would  throw  measurements 
off.  A laser  ranger  would  be  a distinct  possibility,  but  such  devices 
are  expensive.  No  matter  what  the  type  of  sensing,  a reliable  knowledge 
of  height  above  the  guideway  at  two  different  chord  positions  would  be 
all  that  is  required  for  height  and  pitch  determination.  If  relative  and 
vehicle  rates  and  displacements  are  known,  then  the  guideway  rates  and 
displacements  are  but  an  algebraic  sum  of  these  and  easy  to  generate. 
However,  seme  thought  must  still  be  given  to  means  of  sensing  the  vehicle 
and  guideway  states. 

THE  LINEAR  OPTIMAL  CONTROL  PROBLEM 

The  straightforward  optimal  control  problem  involves  finding  the 
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coefficients  of  the  feedback  matrix  which  yield  collectively  the  lowest 
cost  function  value  when  the  cost  function  is  a function  of  the  states 
and  the  control  inputs.  The  desired  cost  function,  however,  is  not  in 
these  terms.  Instead,  the  items  weighted  are  vertical  acceleration,  gap 
variation,  and  control  deflections.  Some  transformation  must  be  made  to 
rewrite  the  performance  index  as  a function  of  actual  state  variables 
and  control  inputs.  Letting  x be  the  vector  of  variables  to  be  weighted 
x be  the  actual  state  variables,  then  the  following  linear  trans- 


formation may  be  made: 


x - M 


which  says  simply  that  the  desired  variables  are  merely  linear  combinations 
of  the  state  variables. 


Since  x is  a vector  of  six  variables  and  x has  seven  and  u has  two, 
M is  a six  by  nine  matrix  of  constants.  Letting  Z be  a diagonal  matrix 
of  weighting  factors  as  previously  used,  we  develop  the  following  cost 
function: 

to 

J»|j  xTZxdt 
*0 

j-U"  u*  *1  mt)  z (m  $)  dt 

to 

J - | f [xT  uT]  M^Z-M  [*]  dt 


The  term  (Ma*Z*M)  Is  the  matrix  inner  product  of  the  transpose  of  M, 
M Itself,  \ Z.  This  product  gives  a matrix  with  off-diagonal  terms 
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and  the  resultant  cost  function  expressed  in  standard  notation  is 


for  a state  system  appearing  as  4 = Ax  + Bu.  However,  it  is  possible 
to  rewrite  this  problem  (Bryson  and  Ho;  1969)  so  that  the  cross  terms 
of  x and  u which  appear  in  the  submatrix  N above  become  zero.  We  first 
rewrite  the  system  equation  as  4 = (A-BP_1NT)  x + Bu,  and  then  the  cost 
function  becomes  _ 


The  submatrices  off  the  main  diagonal  are  now  all  zero  and  the  form  is 
along  more  standard  lines.  For  this  formulation,  we  find  that 

KT  - P"1  (NT  + btr) 

is  the  feedback  matrix  where  R is  the  solution  of  the  matrix  Ricatti 
equation  cited  in  numerous  references  (Schultz  and  Melsa;  1967). 

The  procedure  to  be  followed  starts  with  determining  the  linear 
transformation  matrix  M and  then  selecting  the  values  of  Z,  the  weighting 
factors.  Then  the  matrix  multiplication  must  be  performed  and  the  sub- 
matrices Z,  N 'and  P.  determined.  The  modified  A matrix  A = A - BP~^NT  is 
calculated  as  well  as  the  modified  Q matrix  Q = Q - NP-1NT.  These  are  sub- 
stituted Into  the  matrix  Ricatti  equation  which  now  appears  as 

ATR  + RX  - RBP_1BTR  + Q = 0 

This  particular  form  is  used  since  we  are  interested  in  a terminal  time 
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at  infinity  where  R is  zero.  Several  methods  of  solving  for  R are  available, 


but  the  most  convenient  way  is  to  integrate  backward  in  time  from  the  terminal 


condition  R (®)  = 0,  until  R is  constant  to  within  some  limit.  The  optimal 


feedback  matrix  K is  found  and  we  rewrite  the  system  as  a white-noise- 


disturbed  regulator: 


x = (A  - BK  ) x + Cu 


This  description  of  the  system  is  that  of  which  we  will  analyze  the  ride 


quality  according  to  the  pattern  of  the  previous  chapters. 


This  entire  process  was  programmed  in  AFL  because  of  that  language's 


ability  to  handle  matrices  with  ease.  The  reader  will  find  in  Appendix  E 


the  source  codes  for  the  entire  routine  as  it  was  programmed  for  this  prob- 


lem, along  with  the  appropriate  logic  flow  diagram  of  the  overall  solution 


technique.  The  entire  routine  required  several  cost  function  variations 


before  an  optimal  control  form  was  discovered,  but  the  cost  of  solving  this 


problem  was  less  than  that  for  the  passive  suspension  owing  to  the  greater 


efficiency  of  this  technique.  The  initial  values  of  the  weighting  factors 


was  found  by  exactly  the  same  means  as  that  described  in  the  previous  chapter: 


the  inverse  square  magnitude  of  the  maximum  permissible  value. 


OPTIMAL  CONTROLLER  AND  VARIATIONS 


After  several  trials,  the  set  of  weighting  parameters  which  produced 


the  best  results  were  the  inverse  magnitudes  squared  of  the  following  maximum 


permissible  values 


— n m ~ 


"te  Gap  - °'33  ln- 


- 


^LE  Gap  " 0,33  in* 


It  is  interesting  to  note  that  there  was  no  weighting  necessary  upon  the 


control  surfaces  or  the  pitching  accelerations.  The  control  deflection 


terms  were  free  because  the  gap  control  had  an  influence  upon  how  large 


the  control  deflections  were  and  would  not  permit  large  deflections  which 


result  in  large  gap  variations.  The  controller  would  not  permit  high 


pitching  accelerations  because  over  the  ninety  feet  from  c.g.  to  the  tail 


these  would  necessarily  mean  substantial  gap  variations. 


In  Table  I are  listed  the  values  of  the  feedback  gain  matrix.  None 


of  the  values  are  so  large  as  to  be  unwieldy;  in  fact,  most  are  quite 


small.  The  controller's  response  to  a guideway  pitch  change  of  + 10  is 


to  add  in  a A6  of  + 11.6°  and  a A6  of  - 28.3°  (-14.8  inches  at  the 

r.r  ' 


winglet  tip).  The  pitch  control  to  such  a situation  is  moderate  and 


nearly  exactly  tracking  the  guideway.  The  heave  control  seems  to  be 


leading  the  approaching  grade  by  beginning  to  lift  the  vehicle  earlier  and 


spreading  out  the  upward  acceleration  over  a longer  period  of  time.  On 


the  other  hand,  if  the  vehicle  were  to  have  a one  foot  per  second  rise 


rate,  the  controller  would  move  the  winglet  tips  down  by  0.1  inch  and 


pitch  the  winglets  nose  up  by  0.015  which  corresponds  to  dropping  the 


trailing  edges  by  0.29  inch  and  raising  the  leading  edges  by  0.19  inch 


• 


, 


in  an  effort  to  induce  a negative  moment  to  pitch  the  body  nose  down  and 


check  the  vertical  rise.  This  second  set  of  motions  seems  small  by  com- 


parison to  the  response  to  attitude  changes  and  one  might  guess  that 


ignoring  the  rate  inputs,  while  deteriorating  perf ormance , might  still 


produce  an  acceptable  vehicle.  This  possibility  and  the  results  of  doing 


it  are  discussed  later. 


EIGENVALUES  AND  EIGENVECTORS 


We  now  begin  the  system  evaluation  of  the  actively  controlled  vehicle 


and  proceed  in  the  usual  fashion.  This  configuration  may  serve  as  a 


standard  against  which  any  other  type  of  suspension  may  be  compared,  since 


for  the  given  set  of  weighting  parameters,  this  is  the  optimal  solution. 


Table  I displays  the  eigenvalues  and  eigenvectors  for  the  three  cases. 


Once  again,  as  compared  to  the  rigid  vehicle,  the  actively  controlled 


vehicle  shows  a distinctly  improved  performance.  The  natural  frequencies 


are  considerably  lower  and  the  damping  ratio  for  both  modes  is  around  0.7. 


This  damping  is  somewhat  better  than  the  passive  suspension,  and  much 


greater  them  the  rigid  case. 


The  eigenvectors  show  that  the  low  frequency  — mode  ratio  is  entered 

h 

similarly  to  the  passive  so  that  this  mode  is  primarily  pitch  in  nature 


and  contains  considerably  less  heave  content  them  the  rigid  case.  The 


high  frequency  g-  mode  ratio  is  almost  unentered  from  the  rigid  case. 


Whereas  the  passive  case  high  frequency  mode  is  all  heave,  the  active 


case  contains  a small  amount  of  pitch,  but  this  is  still  so  small  that 


the  high  frequency  mode,  as  in  the  rigid  case,  may  be  considered  nearly 


«.n  heave.  Figure  1 shows  a typical  time  domain  response  of  the  active 


r 


56 


ptgMrr&safcHy- 


TABLE  1 

EIGENVALUES  AND  EIGENVECTORS  OF  ACTIVE  SYSTEM 


Eigenvalues 

-8 

-5.5 

0 

.976  + i 1.014 
.976  - i 1.014 
.408  + i .412 
.408  - i .412 


Eigenvectors 
High  Frequency 


Low  Frequency: 


The  following  are  the  feedback  gains  for  the  controller: 


A6 

A6 

units 

Wl 

w8 

,00474 

.000266 

rad 

fps 

338 

- .0319 

rad 

rad/i 

ACTIVE  CASE  RESPONSE  TO  HOOT • i AND  *•.! 


PLOT  OF  HEIGHT  AND  PITCH i HEIGHT  IS  ST700TH  CURVE 


Figure  1.  Transient  Response  of  Actively  Suspended  Vehicle 
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case,  in  which  the  virtual  separation  of  the  modes  (decoupling  h and  0) 
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can  be  seen. 

ACCELERATION  SPECTRAL  DENSITY 

In  Figures  2 and  3 are  shown  the  acceleration  spectral  densities  of 
the  .40c  and  ,90c  positions  in  the  vehicle.  The  curves  show  a very  smooth 
transition  from  low  to  high  frequencies  and  agreement  with  the  specification 
at  all  frequencies  with  roaa  to  spare.  There  are  no  sharp  peaks  or  valleys, 
only  a very  gentle  rounded  peak  indicating  the  combined  effect  of  the  two 
modes.  When  compared  against  the  corresponding  figures  of  the  other  two 
cases,  the  active  case  proves  to  have  the  lowest  vertical  acceleration  at 
any  given  frequency,  as  one  might  expect.  Once  again  the  reader's  attention 
is  called  to  the  vertical  axis  labelling,  which  is  20  log1Q  I I format. 
The  conversion  to  hertz  is  9.15  rad/r  = 1 Hz. 

RMS  ACCELERATION.  GAP  VARIATION.  AND  CONTROL  DEFLECTION 

The  accelerations  at  the  c.g.  are  greatly  improved  over  the  rigid 
vehicle  (0.121g)  and  somewhat  better  than  the  passive  suspension  (0.0071g) 
having  a rms  value  of  0.0C212g.  This  could  be  considered  a very  smooth 
and  comfortable  ride.  The  rearmost  seat  proved  even  smoother,  being  on 
the  order  of  O.OOlg.  The  seats  in  between  will  have  values  ranging  be- 
tween these  two,  so  that  the  overall  level  of  vertical  accelerations  is 
very  low  indeed. 

The  rms  gap  variation  is  on  the  order  of  one-third  inch,  that  being 
the  criterion  against  which  the  cost  function  was  altered.  There  is  a 
tradeoff  between  gap  variation  and  accelerations,  and  the  cost  function 
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was  adjusted  until  this  level  of  gap  variation  was  met  but  not  exceeded. 
Hence  the  values  of  0.218  Inch  at  the  leading  edge  and  0.327  at  the  trailing 
edge  are  not  so  much  a result  as  they  are  a design  condition  for  the  active 
controller.  These  values  mean  that  the  steady-state  gap  is  more  than  three 
standard  deviations  out  on  the  normal  distribution,  implying  that  contact 
is  rare.  As  with  the  passive  suspension,  there  is  a reversal  from  the 
rigid  case  in  the  relative  size  of  leading  and  trailing  edge  gap  variation. 

The  typical  winglet  motions  are  slightly  larger  for  the  active  con- 
troller; the  rms  A6  is  0.025  rad.  which  corresponds  to  0.75  inch  at  the 

wi 

tip,  and  rms  A6  is  0.00055  rad.  corresponding  to  a trailing  edge  gap 

ws 

variation  of  0.59*+  inch.  These  are  both  significant  motions  with  regard 
to  the  size  of  the  gap,  but  they  are  small  enough  on  the  vehicle  level  to 
imply  that  the  power  required  for  the  winglets  is  likely  to  be  relatively 
small  as  compared  to  total  vehicle  energy  consumption.  The  seemingly  small 
feedback  gains  yield  a fairly  tight,  or  fast-responding,  control  system, 
resulting  in  a smooth  ride  for  the  passenger.  As  with  the  passive  case, 
the  dominant  mode  of  winglet  motion  is  as  a guideway  follower  so  that  both 
gap  and  vehicle  height  remain  nearly  constant.  With  the  active  controller, 
however,  there  seems  an  additional  leading  of  the  guideway  by  the  vehicle 
which  tends  to  spread  the  accelerations  over  a longer  time  period,  helping 
to  smooth  out  the  ride  in  a way  that  the  passive  suspension  could  not  do. 

THE  "PARTIALLY  BLIND"  CONTROLLER 

In  the  process  of  designing  for  cost,  one  might  be  led  to  ask,  "Is 
there  any  way  to  reduce  the  number  of  sensors  and  feedback  circuits  and 
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Partially  Sensing  Active  Suspension  Acceleration  Spectral  Density  at  . 


still  achieve  reasonable  performance?"  The  answer  seems  to  be  yes.  By 
zeroing  elements  of  the  feedback  matrix,  one  effectively  shuts  down  sensors 
and  "blinds"  the  controller  to  those  inputs.  Doing  so  and  performing  a ride 
quality  analysis  yields  interesting  results.  Eliminating  all  sensing  of  the 
guideway  motions  with  respect  to  the  earth  (the  last  three  states)  produces 
no  noticeable  change  in  vehicle  performance.  With  only  vehicle  velocities 
and  relative  displacements,  the  performance  is  almost  exactly  the  same  as 
full  sensing,  the  difference  being  less  than  one  percent  by  any  measure. 

It  seems  that  these  three  of  the  seven  sensors  may  be  eliminated  out  of  hand. 

The  feedback  of  relative  displacements  is  crucial  to  performance  and 
provides  most  of  the  improvement  noted  above.  Eliminating  the  velocity 
feedbacks  does  degrade  performance  noticeably,  but  not  so  far  as  to  be 
intolerable.  With  only  the  relative  height  and  pitch  displacement  feedbacks, 
we  find  that  the  rms  accelerations  increase  to  0.017g  at  the  c.g.  and  0.ce6g 
at  .90c.  While  these  are  not  exactly  considered  a smooth  ride,  they  are 
quite  tolerable.  The  respective  spectral  densities  are  shown  in  Figures  4 
and  5.  Figure  6 shows  the  time  response  of  this  system  to  the  same  set  of 
initial  conditions  as  used  before  and  one  may  see  from  this  plot  that  some 
of  the  pitch  damping  is  gone,  but  that  the  overall  response  is  similar. 

The  indication  is  that  attitude  feedback  is  most  critical,  while  rate  feed- 
back fine  tunes  the  response  and  adds  seme  degree  of  smoothness.  Guideway 
position  and  rate  feedback  is  unnecessary. 

Overall,  the  active  control  does  an  extremely  effective  job  of  smoothing 
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out  the  ride  and  maintaining  a reasonable  gap  clearance.  The  acceleration 
levels  are  somewhat  improved  over  the  passive  suspension,  though  there 
might  be  reasonable  doubt  that  the  improvement  is  worth  the  extra  cost  of 
the  active  system  over  the  passive.  The  active  controller  does,  however, 
yield  accelerations  of  the  order  of  O.OCeg  or  less  (depending  on  position 
inwirip  the  passenger  compartment),  requires  very  small  control  deflections 
and  hence  «man  control  power,  and  allows  us  to  virtually  select  what  levels 
of  gap  variation  the  vehicle  should  have.  The  feedbacks  are  quite  small  and 
therefore  easily  realized.  In  addition,  even  if  the  active  system  is  not  to 
be  realized,  it  does  provide  some  insight  to  the  limits  of  vehicle  performance 
It  seems  that,  as  a result  of  this  analysis,  the  performance  of  the  TRACV  at 
its  best  seems  to  be  very  premising  indeed. 
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CHAPTER  VI 


CONCLUSION 
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In  examining  the  Tracked  Ram  Air  Cushion  Vehicle,  one  finds  himself 
confronted  with  a high  degree  of  heave-pitch  coupling.  This  strong  coupling 
means  that  for  longitudinal  studies,  one  cannot  ignore  one  when  studying  the 
other.  In  addition,  any  control  system  or  suspension  must  be  able  to  avoid 
this  coupling  as  well  as  changing  the  frequency  and  damping  of  the  dominant 
modes.  Using  this  criterion  for  a choice  in  suspension  methods  will  lead 
toward  the  ability  to  positively  control  ride  quality. 

Since  the  TRACV  is  intended  as  high-speed  ground  transportation  for 
passengers,  the  ride  quality  of  the  passenger  compartment  in  terms  of  noise 
and  accelerations  is  of  paramount  importance.  The  vertical  accelerations 
experienced  in  the  passenger  compartment  are  the  subject  of  this  study;  the 
test  subjects  are  the  basic  vehicle  and  two  types  of  vehicle  control  (passive 
element  and  active)  using  two-degrees -of- freedom  winglets.  The  ride  quality 
of  these  three  configurations  over  a surface  having  the  roughness  of  smooth 
highway  is  the  design  criterion  and  the  outcome  is  briefly  described  below. 
CONCLUSIONS 

Using  a linear,  unconstrained  model  of  the  vehicle  dynamics,  the 
following  conclusions  may  be  drawn  about  TRACV  performance: 

1.  The  basic  vehicle  with  no  control  surfaces  is  unacceptable. 

In  almost  all  positions  in  the  passenger  compartment,  the 
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TABLE  1 


SUSPENSIONS  COMPARISON 


Rigid 

Vehicle 


Passive 

Suspensions 


Active 

Suspensions 


RMS  Accel,  at  .4c  (g) 

0.121 

0.0072 

0.00212 

RMS  Accel,  at  ,9c  (g) 

0.0055 

0.0091 

0.00092 

RMS  L.E.  Gap  Variation  (in.) 

14.7 

0.143 

0.218 

RMS  T.E.  Gap  Variation  (in.) 

6.07 

0.224 

0.327 

RMS  A6„  (rad. ) 

" 1 

— 

0.0154 

0.025 

RMS  AS  (rad. ) 

__ 

0.00053 

0.00055 

w8 

High  Frequency  Mode  (Hz) 

1.47 

0.324 

0.224 

Low  Frequency  Mode  (Hz) 

0.962 

0.036 

0.134 

High  Frequency  Damping 

0.226 

0.475 

0.692 

Low  Frequency  Damping 

0.078 

0.337 

0.704 
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acceleration  spectral  density  exceeds  the  UTACV  standard 
in  the  1 - 2 Hz  range.  In  addition,  potentially  destructive 
winglet-guideway  contact  frequently  occurs. 

2.  A two-degree-of- freedom  winglet  suspended  by  springs  and 
dampers  smooths  out  most  of  the  accelerations  of  the  rigid 
vehicle.  Guideway  contact  is  rare,  the  accelerations  meet 
the  specifications,  and  the  modes  of  motion  are  almost  com- 
pletely decoupled  in  addition  to  the  lowering  of  vehicle 
natural  frequencies  and  increasing  damping  by  allowing  the 
winglets  to  become  guideway  followers. 

3.  The  active  control  of  the  two-degree-of- freedom  winglets 
further  improves  the  performance.  The  accelerations  reach 
a level  approaching  the  threshold  of  sensitivity  while  not 
requiring  large  amounts  of  power  to  do  so.  This  configuration 
however,  gives  only  moderate  performance  improvement  over  the 
passive  suspension  (which  is  quite  acceptable)  in  return  for 
substantial  construction  and  maintenance  cost  increases. 

4.  If  the  active  controller  is  to  be  used,  the  sensing  of  the 
guideway  motions  is  irrelevant  to  ride  quality.  Rate  sensing 
while  not  strictly  necessary,  does  give  substantial  improve- 
ment at  levels  where  improvement  is  needed. 

DIRECTIONS  FOR  FURTHER  INVESTIGATION 

One  point  of  this  study  subject  to  question  is  the  guideway  approximation. 
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The  average  of  guideway  pitch  and  height  is  a fine  descriptor  for  wave- 
lengths much  smaller  or  much  larger  than  the  length  of  the  vehicle.  In 
the  transition  region  near  where  the  length  of  the  vehicle  is  about  one- 
half  the  guideway  wavelength,  the  effects  not  accounted  for  in  the  second 
order  linear  approximation  mi./  be  of  importance,  particularly  with  regard 
to  the  effect  on  moment  forces  of  such  a bump  as  it  traverses  the  length 
of  the  vehicle;  these  bear  further  study. 

Different  schemes  for  control  also  seem  to  be  a likely  prospect  for 
study.  Mechanization  of  a two-degree-of-freedom  winglet  may  be  no  small 
task,  and  there  might  be  alternatives.  One  such  possibility  is  pitch 
control  by  way  of  a trailing  edge  flap.  This  possibility  requires  study 
of  the  dynamics  of  a two-segment  bottom  slope.  Another  possibility  is 
that  given  the  pitch  control,  A6  , is  very  small,  allowing  the  winglet 

wa 

only  A6  freedom  and  warping  the  winglet  to  obtain  A6  effects  may  be 

w3 

a distinct  possibility.  This  prospect,  however,  requires  close  study  of 
stress  and  fatigue  factors. 

Finally,  vehicle  dimensions  will  play  a strong  role  in  the  dynamics. 

A vehicle  fifty  feet  long  is  as  likely  as  one  one  hundred  fifty  feet  long. 
Top  surface  lift  will  have  an  effect  on  moment  equations.  The  flat,  sloping 
bottom  surface  is  an  important  part  of  the  present  model  as  well  as  aspect 
ratio  and  weight.  Since  all  of  these  factors  will  significantly  affect  the 
dynamics,  it  is  Important  that  a practicable  vehicle  be  laid  out  before  any 
further  detailed  work  on  longitudinal  ride  quality  is  undertaken. 
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APPENDIX  A 
MJATIONS  OF  MOTION 


The  model  of  the  TRACV  used  in  this  thesis  is  a two-degree-of -freedom 
aerodynamic  mode  with  freedom  in  heave  and  pitch.  It  assumes  that  the 
lateral  and  longitudinal  aspects  may  be  considered  separately,  and  the 
present  concern  is  on  longitudinal  motions  only.  Other  assumptions  in- 
clude constant  forward  velocity,  constant  mass,  and  fixed  center  of  gravity. 
The  lift  and  moment  may  be  described  by  Taylor  series: 

L = L + Ah  + Aq  + — Ah  + §£  A0  + (higher-order  terms) 

° ah  3q  ah  30 


M = M + Ah  + |^  Aq  + — Ah  + !£  A0  + (higher-order  terms) 

° ah  3q  a&  30 

This  expression  of  the  lift  and  motion  reflects  the  basic  vehicle  dynamics. 
Assuming  small  perturbations  of  the  motion  quantities  only  around  some 
steady-state,  the  higher  order  terms  may  be  neglected.  In  steady-state, 

M *»  0,  L = mg,  and  there  are  no  perturbations,  so  we  may  write: 


mg  = L = L 


0 = M = M 


The  sum  of  all  the  steady-state  terms  is  zero  and  they  may  be  dropped 
along  with  the  higher-order  terms.  In  addition,  if  the  winglets  are 
allowed  to  move,  then  the  lift  and  moment  have  derivatives  with  respect 
to  the  control  motions.  In  final  form  the  equations  are: 


AL  » ^ A&  + jgf  Aq  + & Ah  + §&  A0  + 

3ft  3q  dh  30  36w1 


(A-l) 
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au  - an  aI  4.  a**  Ar  . m . sm 

AM  = Ah  + g-  Aq  + — Ah  + gg-  A0  + 

oh  oh  Wj 


(A -2) 


It  should  be  noted  that  the  motion  variables  above  are  all  of  the  vehicle 
relative  to  the  guideway. 

Having  expressed  the  change  in  lift  and  moment,  we  may  now  write 
relations  concerning  vehicle  accelerations: 

mE  = AL 

v 


mch  =rp  U3S  ACt 
v 2 r L 


- — — h - ACt 

5 p U8  S V L 


h = ACt 
v L 


(A-3) 


CCL 

The  term  has  the  units  sec3  so  that  one  convenient  way  to  obsorb  it 

is  to  change  to  non-dimensional  time.  Such  a conversion  is  convenient  in 
that  if,  during  the  scaling  process,  the  term  ^ is  kept  constant  (A  defined 
below),  the  dynamics  will  be  unchanged  except  for  frequency.  We  then  write: 


seconds 


t 

T “A 


(A-4) 


In  making  this  conversion,  one  must  be  careful  that  any  time  related  terms 

•• 

be  adjusted.  All  the  time  derivatives  must  be  rewritten  as  is  h below: 


""  lit  mtgjiaaadugiaah 


- ,»a  _ *a  _ 

B - 55  “ " dOT  h * a*  h 

equation  (A-3)  can  be  rewritten: 


L£  (i  i)  - A*  (i  B*  ) - S'  - AC 


8 Aa 


Remembering  from  Chapter  II  the  relationships 


t c * 

h = uh 


*-§< 


we  are 


able  to  express  equation  (A-l)  in  the  following  forms: 


c,-  45  + ± CL  40'  + CTi.  45  ♦ Cu  40 


+ °L  «»,  + % 

w,  w8 


(A-5) 


The  compact  notation  C_  means  -7-  and  so  on. 

Lh  dh 

Following  a similar  argument  yields  a moment  equation: 

e — C Ah'  + ^ Cu  A0'  + CL.  Ah  + C A0 
v UA  MT  1 M ”0 


+ \ 


AS  + C„  A6„  AS 


'1 


w8  w8 


k a 

In  the  development  of  (a-6)  however,  the  term  (-^)  comes  up,  where  k^ 


ca  9CM 


ia  the  vehicle  radius  of  gyration.  We  define  CM_  to  be  ^ g and 


(A-6) 


so  forth  to  absorb  it.  In  the  present  model,  k^  was  taken  to  be  the 
value  of  the  gyration  radius  of  a uniform  rigid  bar  of  vehicle  proportions, 
Adding  the  two  relationships  h - b and  ©'=  0 to  (A-5)  and  (A-6) 
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gives  us  our  fourth  order  pitch-heave  model. 

The  steady-state  lift  coefficient  is  CL  ■ .455,  a value  found  to 
be  typical  of  those  found  in  towed  model  tests  at  the  Princeton  University 
Dynamic  Model  Track.  With  an  assumed  forward  velocity  of  300  fps,  the 
vehicle  weight  for  equilibrium  is  110, 000  lbs.,  agreeing  fairly  closely 
with  the  Mitre  vehicle  (Fraize  and  Barrows,  1973).  Also  the  present 
model  implies  A = 1.456  sec.  It  should  be  noted  that  this  theory  assumes 
no  upper  surface  contributions  to  either  lift  or  moment.  In  actuality, 
the  shape  described  adds  a ACL  = .25,  which  will  slightly  alter  moment 
characteristics  and  allow  a heavier  vehicle. 
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APrtiHDIX  B 

STABILITY  DERIVATIVES 

Presented  below  are  the  equations  for  the  attitude  stability  de- 
rivatives as  they  were  theoretically  developed  in  (Curtiss  and  Putman; 
lSyTT).  The  rate  derivatives  are  the  result  of  a quasi-static  theory 
developed  by  Professor  H.  C.  Curtiss,  Jr.  of  Princeton  University,  and 
are  based  on  the  under-vehicle  pressure  which  assumes  a level  pitch 
attitude  (0  * 0)  and  incompressible  flow.  The  effect  of  vehicle  rates 
is  as  an  increment  to  the  lift  and  moment  coefficients. 
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3h 
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5 8 3CL  . 2c  aCL 
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- ,1  r ?V  t'fil 
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^s-  (J  - *-±i)  ^ - - (i  +tn 

8h  5 a&  As  53  1+5 
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aq 


- . , acT 

/-  cr  + lx  L . 

(lto--I->3T  + 


M°  [2  (1  + 7 - x )(S  + tn  — ^_) 

A*  or3  a 1 + a 


- -c|  ([1  + 5]3  -1)  +tn  — 
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The  various  partial  derivatives  with  respect  to  the  ratio  parameters  may  be 
found  in  the  cited  reference  and  were  computed  by  the  FORTRAN  program  at 
the  end  of  this  Appendix. 

To  find  the  derivatives  of  lift  and  moment  derivatives  with  respect 

to  the  control  surface  deflections,  we  express  them  as  effects  on  the 

dimensionless  parameters.  We  find  that  6 affects  only  r , and  6 affects 

Wj  o’  w8 

only  rx.  Figure  B-l  illustrates  a perturbed  winglet  proximate  to  a guide- 
^ way  lip.  Examination  of  this  geometry  shows  that,  assuming  a guideway  lip 
cant  of  forty- five  degrees  and  ancminally  level  winglet,  on  can  write 


A60  = .707  W,  i A5Wi 

and,  incorporating  the  definition  of  rQ,  one  can  express  the.  lift  coefficient 
and  moment  coefficient  derivatives  as  follows: 
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II 
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B5704.  DOT7R5IN,TI?{E  = 59,  PAGES=260:,  C=0 
UNCIIED  JOB  OPTION.  C=0  IS  INVALID 


$JCE 

HNI!IG**  MS 


3 El S PROGRAM  CALCULATES  THE  CHA 5ACIERISTICS  OF  A PAM-AI 3 
CUSHION  VEHICLE  BY  DETE3U IllIA'G  THE  CHARACTERISTIC  EQUATION 
COEFFICIENTS,  THE  SOOTS,  AND  THE  SI ABILITY.  DERIVI7IV33. 

THESE  ARE  CALAULATSD  OS  PCU3  INPUT  CARDS.  ' 

I EE  FIRST  CARD  SHOULD  EE  INITIAL  PITCH  ISfliE  AMD  INITIAL 
2 SUE  ZEECL  IS  2?  10. 5 FORMAT.  THE , SECOND-  CARD  SHOULD,  RE  CHORD  AND 
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E CARD  SHOULD  B 


WIETK  IN  2F10.5  FORMAT.  THE  THI 
ALPHA  EAR. 

AIL  ANGLES. AEE  IN  RADIANS,  El 
IT  IS  PCSSI3LE"T0  ITERATE  MCR 
BY  PIACING  VALUES  FOR  FLAG, LIMIT 
CASE  IN- II  ,F10.5,F10.5  FORMAT.  .. 

IF  NC  ITERATION  IS  DESIREE . CO 
"1"  = ALPHA  BAR 

”2 11  = BOTTOM  SLCPE  

f;-  CHORD  --  W 

*-I SITI&L  PITCH  ANGLE'-. 

sub.  .zero.  . 

"6"  = WIDTH 

INITIAL  VALUE  FOR  THE  ITERATION  IS  TH 
FIRST  THREE  CARDS.  FINAL  VALUE  IS  3E 
ITERATIVE  PROCEDURE  STEPS  ONLY  UPWARD 


THIS  SEGMENT  GIVES  THE  STABILITY  DEEIVA7XV2S  AND 
"PA»  KEANS  THETA  (PITCE  ANGLE) 

"ECTTOM”  MEANS  ALPHA  Sl'B  ZERO 
EXCEPT  IN  "BOTTOM",  ANY  »C"  CHARACTERS  ARE.  ZEROES 
INPUT  DATA  ARE:  CHORD  (FEET) , THETAO  (RADIANS) 

WIDTH  (FEET),  HC,  BOTTOM, 
ALFBAR,  XO. 


IMPLICIT  S EAL*8  (\-H,0-$) 

EXIEENJL  F 1 , F2 

C0MPLEX*16  ROOT  1 , SOOT 2 , CD  SORT , ROOT 

REAL  LIMIT 

INTEGER  FLAG 

REAL  RC,K, THETAO, BOTTOM 

A ICG  |X)  = CL CG(X) 

AT  AN  (X)  = DAIAN  (X) 

SCB T IX  ) = DS  CBT  ( X ) 

SEAL  15,150) C,W  _ : 

REAL  (5, 150)BOTTCM,  ALFBAR 
ECECAT  (2F1C.5) 

REAL  15, 151)FLAG, LIMIT, STPSIZ 

FCEfcBT  (I1,F10.5,F10.5) _ 

TBFTlC  = 0 - -V  rr--  7 

BTS*2.  ' - V 

DEI C=HTE-2 .916  

AI FEAR -C/HTE*  (THETAO +BOTT CM) 
BC=2*EEL0/ (W^(THETA0+ROTTCM) ) 
Fl*-2*C*THETA0/  (W*  (TH  ETAO  + L’OTTOM)  ) 
AIFCFG=ALFEAR 
ECTCFG=9CTT0M 

HCBG*W  . _ il. 

CCFG=C 
SCCE  G=BC 

IHECEG*THETAO  8l 


SOOTC, 300TD 


Etll?tA=SCHI  ( (1 .0+.UFBAR)  **2-1.  C)  DCCT’’ 

F21=-2 .0*£1A»DA/SCHT  (ELA*  EA  **2  *1  .0)  |Jf\J 

1 + AI*MELAMDA)/SCHT(ELAHDA**2  + 1 .0)  . ' ..■f.fSfir 

2 ♦ AICG  (SLARCA+SC?? (ELAHDA**2+1 .0)  ) 
F3L=FLAHDA-2.0*A?Ay  (ELAKDA)  **105  (ZL1RDA+SQ5Z  ( 

1 /SCFT  (ELA!!DA**2*1 .0) 

F2H=£L AHDA-ATA3  (ELAMDA)  -BKTGEI  (Q.0EC0,£LA 

F3a=Fi*aoA+scfi?  (ela:ida**2*i  .0)  *.  5 

1 -.S’*  A LOG  ( ELAMDA  + SQET  (ELA EDA**2+ 1 . C)  ) 

2 -5NSGFL(0.0D00, ELAMDA,  F2  ,.0G1D00) 

£12=2. C+F2I/ALFEA3 

E12=  (F3L-2.0*F2L)/ALFBAE**2  . .1 .. 

ER2=2.C*F2  S/ALFBAH**2 
EK2  = (F3M-2.0  + F2M) /ALFBA3**3  ;•  ‘ 

D1EA=  (1.0  + ALFBAE)  *(  ( 1 . 0+ A 1FBA3)  .**2-  1 . 0)  .** 

E 2 1 1 = (1./SCaT(ELA!!lDA*2LAMEA+1.)*  (-2.0)  +EI 
1*  (SIAaLA*ELARDA  + 1 .)  **  (-1 . 5)  ) *DLDA 
F 2 1 2=  ( 1 . 0/..CELA M D A*EL  A BD A*  1 . )-**  .3*1 . 0/  ( EIA 
1 ATSK  (fLABDA):*  (-1.0)*  (ELA  KDA**2+1 . 0)  **  (-1 
F 21 3 = ( (I.C+ELAMBA*  (E1ABDA**2  + 1 .(£)**  (-  .5) 

1. I . 5)  ) ).*DLD  A - -li 

DF2l=F2L1+F2L2+F2L3 
F 21 1=1 .0*DLDA 

F3l2=-2. 0/  (ELA!iBA**2 +1 .0)  * D I D A 

F3I2=(1.0/  (ELAHDA**2+1.Q)  **  (.5)*  (1.0+ELAa 
1 **  (-.5)  ) / (BLAMDA+  (5L  AIiDA**2+T.  0)  **  (.5)  ) 

2+AICG  (FLAMEA+  (2LAMDA**2/H  .0)  **  (.5)  ) * (-1.3 
3**  (-1 .5) ) ) *DLDA 
EF3I  =F3L1+F3L2+F3L3 

F2Ml=CIDA  . . . 

F2K2=- 1 . 0/ (EXAM BA**  2+ 1 .0) *DLDA 
F2E2= (-ELA  EDA*AXAS (ELAHDA ) / (ELAM DA**2+ 1 . C 
EF  2K=F  2M 1 +F2M2+F2M3 

F3l11=(.5*  ( FLAME  A**2  + 1 . 0)  **  ( . 5)  +.  5* FLAME A* 

1 )*EIDA 

F3ri2  = -.5*  (1.  Q-*- ELAMDA*  (ELA KDA**2+1 . 0)  **  (-. 

i)  **  (.5) ) *elea 

F3Sl3=-  (ELA  MDA*ALOG  (ELAMDA + (ELAMDA**2+1 . 0) 

1 *EIIA  _ . . -,:aV; 

Df 3£=F2al+F3H2+F3M3 

EL  2 E ?.=  2. 0*  (1 . 0/ALF3AR*EF2  I - 1 . 0/ALF DAH**2 
Dl3Ei=1. 0/ ALFBAS**2* (DF3L-2 . 0* DF2L ) + (-2 .0 
1 2 . C*F2L)  - - WM& 

DM2EA=2.0*  (1 .0/ALF3A3**2*  EF2K-2. C/AIF3AR* 

DS3  E£=1 . 0/ ALFBAS**  3*  (DF3fi -2 . C*DF2K)  - • 
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AFTENDIX  C 

FINITE  PAD  LENGTH  EFFECT 

As  on  air  cushion  vehicle  travels  over  the  ground.  It  encounters  a 
series  of  surface  irregularities  which  alter  the  conditions  of  the  air 
cushion.  Since  the  irregularities  occur  in  all  wavelengths,  it  is  likely 
that  sane  of  the  irregularities  will  be  of  much  shorter  wavelength  than 
the  length  of  the  air  cushion.  When  this  is  the  case,  the  shape  of  the 
ground  with  respect  to  the  vehicle  cannot  be  described  by  a single  number, 
as  is  the  case  with  a point  contact  suspension,  but  is  instead  some 
function  of  the  pad  length,  or  length  of  the  air  cushion.  For  instance, 
to  a point  contact  suspension  a stairstep  is  a step  change  in  surface 
height,  but  to  a pad  of  finite  length  it  is  an  obstacle  to  be  climbed, 
for  a step  change  in  surface  height  implies  that  all  points  under  the 
pad  change  height  simultaneously,  which  is  not  the  present  case. 

The  model  in  this  paper  uses  attitude  and  rate  terms  to  describe 
the  effects  of  vehicle  motions  on  the  forces  on  the  body.  Because  of 
this  fact,  the  typical  surface  spectrum  will  cause  the  vehicle  to  have 
a constant  acceleration  spectral  density  high  frequency  asymptote  when 
constant  guideway  descriptors  are  used.  This  implies  infinite  rms 
accelerations  which  is  obviously  wrong,  and  experience  shows  that  very 
high  frequencies  are  usually  severely  attenuated.  The  discrepancy  is 
caused  by  the  model  chosen  which  implies  that  the  relative  height  and 
pitch  of  the  vehicle  over  the  guideway  may  each  be  described  by  a single 
value.  In  actuality,  where  irregularity  wavelengths  on  the  order  of 
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and  shorter  than  vehicle  length  are  present,  the  relative  height  and  pitch 
are  necessarily  functions  of  distance  along  the  vehicle. 


Including  such  functions  in  the  model  make  any  kind  of  analysis  ex- 
tremely tedious  if  not  impossible.  Simplification  is  possible,  however. 
Suppose  we  assume  two-dimensional  irregularities  such  that  the  surface  is 
constant  across  the  track  width  at  any  point  and  that  the  average  height 
and  pitch  of  the  guideway  portion  under  the  vehicle  are  the  most  significant 
descriptors.  It  then  becomes  possible  to  describe  the  effect  of  vehicle 
size  on  small  disturbances  as  an  attenuation  by  1 (Riblich,  Captain, 
and  Richardson,  1967). 

For  ease  of  description,  we  assume  at  tQ  = 0 the  center  of  the  vehicle 
is  at  a crest  of  the  guideway  and  that  the  guideway  itself  is  expressed 
as  a sum  of  sinusoids  over  a continuous  frequency  domain.  Using  X to 
represent  guideway  wavelength,  the  height  of  the  guideway  with  respect 
to  sane  mean  datum  plane  is 

Sg-|Bg(x)|  cos 


where  x is  distance  in  front  of  the  vehicle  center.  If  we  normalize 
guideway  amplitude  to  one,  then  the  average  height  of  the  guideway 
under  the  vehicle  is 


and  by  trigonometric  identity 


z _ / X ttL\  2rrUt 

Vvg  ^ ,la  r> cos  — 


(C-l) 


The  term  in  parentheses  in  equation  (C-l)  is  the 


function  for 


5 " “•  This  function  is  plotted  in  Figure  C-l  which  demonstrates  that 
this  function  goes  to  zero  for  very  short  wavelengths.  This  result  is 
not  surprising  since  experience  tells  us  that  disturbances  with  a wave- 
length of  one  inch  will  be  of  no  effect  on  a large,  150  foot  long  vehicle. 

The  rate  of  guideway  height  change,  h,  is  found  in  a similar  fashion. 
Normalized  height  was  expressed  by 

h = cos  <S2£  + ) 

g X X ' 


so  its  time  derivative  is 


-2ttU 


The  average  value  is 
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This  last  expression  is  the  time  derivative  of  (C-l). 

It  turns  out  that  the  effect  of  averaging  guideway  pitch  is  similar. 
We  proceed  in  the  following  way 

dh 

t“ee-s£ 


tan  0 


g dx 
dt 


eg  - g 

For  the  guideway  model  we  have  chosen,  as  with  the  vehicle,  the  small  - 
angle  assumption  is  valid.  This  being  the  case,  one  can  say  tan  9 - 0 
and  as  a result 

A 

, . h 

g U 

It  then  follows  that 


f-"  <&  ■«“  T> 


(C-3) 


and,  in  a similar  manner 


e - iali!  (i.  sin  TSl)  cos  (0-4) 

*Avg  "L  1 x 

It  is  apparent  that  equations  (C-2),  (C-3)»  and  (c-4)  are  related  to 
(C-l)  by  constants  and  time  derivatives. 

If,  as  we  have  seen,  the  effect  of  using  average  guideway  parameters 
is  to  add  attenuation  to  the  guideway  effects  as  a function  of  frequency, 
then  this  attenuation  must  be  Included  in  the  overall  model.  It  is  possible 
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to  approximate  the  function  1 ^ I with  a linear  second  order  filter* 

Using  a linear  filter  gives  us  the  advantage  of  freedom  in  incorporating 
it  anywhere  into  the  model.  A second  order  filter  is  necessary  to  be 
assured  that  the  effect  of  height  and  pitch  variation  in  the  guideway 
will  diminish  with  increasing  frequency.  Figure  C-l  shows  the  actual 
averaging  attenuation  function  superimposed  on  the  transfer  function 

" (s  + 8)(s  + 12) 


Where  s is  the  Laplace  operator  in  non-dimens ional  time.  This  illustration 
depicts  the  close  agreement  between  the  linear  approximation  and  actual 
function. 

As  previously  stated,  the  height  and  pitch  spectra  of  a typical 
guideway  are  given  as 


The  apparent  guideway  to  the  vehicle  is  further  attenuated  by  the  size 
of  the  vehicle.  One  may  state  then  that  the  apparent  guideway  spectra 


as  the  vehicle  sees  them  is  given  by 

„ / x 96\^U 

g'  ' s (s  + 8)(s  + 12) 


(C-5) 


96 

*«(,)  " (•  ♦ 8)(»D+  12)  (c_6) 

Equations  (C-5)  and  (C-6)  represent  the  guideway  spectra  as  they 
appear  to  the  vehicle.  This  must  be  part  of  the  overall  model  so  that  the 
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controller  and  dynamics  analysis  may  account  for  it.  These  two  equations 


1- 

represent  a close  linear  approximation  to  a function  which  is  the  result 
of  assuming  that  the  parameters  of  primary  importance  are  the  averages  of 

K! ' 

' height  and  pitch  under  the  vehicle  at  a given  point  in  time.  These  can 


be  considered  as  a filter  placed  on  white  noise  so  that  the  source  of 


disturbance  for  control  theory  purposes  is  white  noise  although  the  guide- 


way itself  is  not,  allowing  the  controller  to  account  for  non-uniform 


guideway  spectra  without  complicating  the  analysis . 
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STATE  MATRIX  FORM 

i 

The  most  convenient  way  to  express  the  equations  of  motion  is  in  the 
state  matrix  form.  From  Appendix  A we  draw  (A- 5)  and  (A-6)  as  the  first 
two  equations.  This  pair  represents  a fourth-order  coupled  dynamic  system 
solved  explicitly  for  vertical  and  pitching  accelerations.  This  pair  has 
a form  which  is  not  adaptable  to  the  state  variable  form  directly,  so  we 
require  the  pair: 


~ (hj  = S' 

dr  r r 


which  defines  the  vehicle-guideway  relationships. 

Appendix  C develops  the  apparent  spectra  of  guideway  height  and  pitch 
variations.  If  7] (r ) represents  the  white  noise  disturbance  of  unity  spectral 
density  (using  the  non-dimensional  time),  then  the  apparent  guideway  differ- 
ential equations  are: 

V + 20  V + 96  eg = 96  T|(t) 


where  the  prime  notation  denotes  the  differential  with  respect  to  t,  non- 
dimensional  time. 

The  guideway  relationships  sum  up  to  a third  order  set,  and  the  vehicle 


variable  form  and  the  state  variables  are:  Ah  ' , A9  ' , Ah  , A0  , 

r r r r 

Ah  ' , A0  A0  > the  control  variables  are : AS  , AS  . These  equations 
g g g wi  wa 

are  expressed  in  numerical  form  in  Table  D-l.  This  controllable  set  of 


equations  for  the  actively  suspended  vehicle  may  be  considered  the 
rigid  case  if  all  elements  in  the  B matrix  are  set  to  zero. 

Table  D-2  shows  the  generalized  form  for  the  passive  suspension 
and  lists  the  substitutions  which  are  not  part  of  the  active  model. 


This  Table  shows  the  modifications  to  the  equations  as  discussed  in 
Chapter  IV,  with  A6  and  A6  now  appearing  as  state  variables  and 

W1  W3 

the  B matrix  eliminated.  From  D-2  one  may  see  the  generalized  form  of 
the  contents  of  D-l  by  striking  the  last  two  rows  and  adding  the  second 
term  which  appears  as 


TABLE  D-l 


NUMERIC  MATRIX  REPRESENTATION  OF  TRACV  STATE  EQUATIONS 
(ACTIVE  CONTROL  FORM) 
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MATRIX  REPRESENTATION  OF  TRACV  STATE  EQUATIONS  (PASSIVE  SUSPENSION 
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APPENDIX  E 


OPTIMIZATION  SOLUTION  FORMAT 

The  optimization  of  the  passive  and  active  suspensions  was  Implemented 
on  the  digital  computer  using  the  interactive  language  known  as  APL.  The 
two  techniques  varied  somewhat  from  one  another.  Figure  E-l  shows  the 
logic  flow  for  the  passive  case.  Figure  E-2  depicts  the  flow  for  the  active 
case.  Following  these  is  a set  of  listings  of  the  source  codes  used.  The  * 

first  listing  is  in  FORTRAN  and  solves  the  mean  squared  error  integral. 

The  remainder  are  in  APL. 

The  APL  Programs  are  designed  to  find  the  numerators  of  the  various 
transfer  functions  of  interest  and  convert  them  to  rms  values  or  sets  of 
poles  and  zeros  for  the  spectral  density  routine,  which  is  a modified  Bode 
plot  routine. 

COSTFN  - performs  the  linear  transformation  which  sets  up  the  per- 

formance index  in  terms  of  state  variables  and  inputs  and 
converts  to  standard  notation  as  discussed  in  Chaper  V. 

FRAME  - converts  the  A matrix  to  a vector  containing  the  coefficients 

of  the  characteristic  polynomial. 

G5AP  - finds  the  numerator  coefficients  for  the  L.E.  and  T.E.  gap 

♦transfer  functions. 

MAKE  - generates  the  A matrix  for  the  passive  suspension  given  the 

values  for  the  spring  and  damper  constants. 

19 
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MEANSQ  - calculates  the  mean  squared  error  integral  for  an  eighth 
order  transfer  function  for  the  passive  solution.  Answer 
given  is  square  root  of  mean  squared  value. 

NUMB  - (and  NUM  for  Passive  System)  finds  the  numerator  polynomial 

for  any  of  the  state  variables. 

OPT  - performs  the  four-variable  gradient  search  for  the  spring 

ftnH  damper  constants  according  to  the  performance  index. 

QUAD  - finds  the  rms  values  of  L.E.  gap  and  T.E.  gap  as  well  as 

rms  control  deflections  for  the  performance  index  use  in 
OPT. 

PIC  - solves  the  matrix-Ricatti  equation  by  reverse-time  integration. 

ROOTS  - finds  the  roots  of  a polynomial  by  calling  the  Bairstow  Method. 

SEAT  - finds  the  rms  vertical  acceleration  at  a given  position  in  the 

vehicle,  given  as  the  decimal  portion  of  the  distance  from 
front  to  rear. 

SET  - arranges  the  numerator  polynomials  for  gap  and  control  varia- 

tions and  calls  MEANSQ. 

VECTOR  - calculates  eigenvect'  s for  complex  eigenvalues. 

WING  - finds  the  rms  contr  ‘flections. 

Particular  attention  must  be  to  polynomials  having  a zero  constant 

0 

term  which  sometimes  do  not  conver  in  the  Bairstow  routine.  When  this  occurs, 

a convenient  method  is  to  drop  it  and  insert  the  root  at  the  origin  later. 
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Determine  Objective  Function 
Weighting  Parameters 


v ^ - $JOB  B1002. DOTTR  AIN,  TIM  S*20,  PAGES* 

**R ARMING**  MIS-PUNCHBD  JOB  OPTION.  BIN=G  IS  INVALID 

£********************♦************************************************** 

C**  '** 

c**  THIS  PROGRAM  CALCULATES  THE  MEAN  SQUARED  ERROR  INTEGRAL 

C**  POR  A TRANSFER  FUNCTION  WHICH  IS  OF  N ORDER  IN  THE  DENOMINATOR  AND** 
C**  f N — 1 ) ORDER  IN  THE  NUMERATOR.  THIS  INTEGRAL  HAS  BEEN  SOLVED  IN  ** 

C**  CLOSED  FORM  AND  IS  PRESENTED  IN  SEVERAL  TEXTS  ON  LINEAR  CONTROL  ■ * 

C**  DESIGN.  THIS  PARTICULAR  SET  OF  EQUATIONS  HAS  TAKEN  LITERALLY  FROM** 
C**  NEWTON, G .C. ,L. A. GOULD , AND  J. F.  KAISER,  ANALYTICAL  DESIGN  OF  LINEAR* 
C**  FEEDBACK  CONTROLS,  1957,  PP.  371-377.  * 

C**  THE  INPUT  FOR  THIS  PROGRAM  TAKES  THE-FOLLOWING  ORDERtr  7 H.  .** 

- c**  THE  : FIRST  CARD  CONTAINS  THE  NUMBER  OF  TRIALS  TO  BE  RUN  (12  FORMAT)  ** 

ts#  c**  THE  SECOND  CARD  HAS  THE  TYPE  OF  SYSTEM  (12, 7X)  AND  A MEMO  ** 

=Ct;  i.  c**  STARTING  AT  COLUMN  10  * >">  CONTAINING  SIXTY  ALPHABETIC  CHARACTF.RS  **j 

C**  THE  NEXT  CARD  IS  THE  'r0R  COEFFICIENTS  IN  DESCENDING  ORDER  ** 

C**  IN  (N)  FI  0.4  FORMAT  A ' AFTER  IS  THE  POLES  IN  (N*1)F10.4  ** 

C**  UP  TO  EIGHT  DA 1A  - LL  FIT  ONTO -A  CATD;  IF  THERE- ARE-  - ** 

^ C**  MORE,  THEN  THE  DATA  "Hu,  LIT  BETWEEN  TWO  CARDS  AS  EVENLY  AS  ** 

C**  POSSIBLE.  IN  THE  EV3UT  C A*  ODD  NUMBER,  THE  FIRST  CARD  GETSONE  ** 
i'V- C**  MORE  THAN  THE  SECOND.  - - 

c**  ANY  GAINS  TO  BE  INCLUDED  MUST  FIRST  BE  DISTRIBUTED  TO  ALL  THE  ** 

. C**  COEFFICIENTS.  ** 

c**  FORTRAN  PROGRAMMING  BY  B.  A.  LUHRS^-  - 

cit**^*************************************?********?*******?********** 

1 IMPLICIT  RSAL*8  (A-H,M,0-Z)  - i.r  2 , - •:!  J 

2 DIMENSION  MEMO  (15) 

3 PRINT,  • THE  RMS  VALOE  IS  GIVEN  FOR  EACH  INPUT.  « 

4 . - READ  (5,201)  ICOUHT  . - ~ -v.c 

; 5 : _ ; DO  999  I J=  1,1  COUNT  ' Y-' Vr:  ~ - -rT  "Y  *v " • "><1  C -&&&£&& 

: 6 READ  (5,100)  I0SDEB,(  ME  HO  (WIN)  ,NIN=1,15)  ' 

7 100  POR  MAT  (I2,7X,  15  A4)  r 

8 104  FORMAT  (4F10.  4) 

9 105  FORMAT  (5F1 0.4) 

10  106  FORMAT  (6F1 0.4)  - --fc 

11  107  FORMAT  (7F10. 4)  - \ . .YYfi'  • . • J 

12  108  FORMAT  (8 PI 0.4)  - 7 

13  200  FORMAT  (E12. 4)  - - • ----- 

14  201  FORMAT  (12) 

15  301  FORMAT  (IX,***  NOTE:  *,15A4) 

16  IF  (IORDER.EQ.7)  GO  TO  7 ~r  : • *•  - c-.r  ^ ■■■' 

17  IF  (IOBDER.  EQ.  8)  GO  TO  8 ' " V >;77:r-7 

18  IF  (IORDER.EQ.  9)  GO  TO  9 ; - Y 

19  6 CONTINUE 

C»** ************** ************************************************* ***** 
^ *« 
c** 

C**  THIS  SEGMENT  PERFORMS  THE  SIXTH  ORDER  OPERATION.  •-  *-• 

C** 

£***«*  ************.*♦*******•*********•**•***♦***••*******  *************** 

20  BEAD  (5,1  06)  C5  ,C4  ,C3,C2,C1,C0  .v  -•  >i« 

21  BEAD (5 ,107)  D6,D5,D4,D3,D2,D1,D0 

22  M1*-D0*D1*D5  *D0*D3**2  *D1  **2*D4-D1  *D2  *P3 

23  H2*D0*D3*D5*D1  **2*D6-D1*D2*  D5 

24  M3* DO*  D5**2*D1*D3*D6- D1*D4*D5 

25  M4* (D2*H3-D4*W2*D6*M  1 ) /DO  > r \ .yi 

26  M5*(D2*M4-D4*H3*D6*H2)/D0  t*: 

27  H0*(D4*H1-D2*H2*D0*M3)/D6 

28  DEL6*  (D1*H5-D3*H4*D5*M3)  *D0  _ 

29  BB?G3l»(C5**  2*H0  ♦ (C4**2-2*C  3*C5)  *M  1*  (C  3**2-2*C2*C4  *2*C1*C5|  *M2* 


1 (C2**2-2*C1*C3  *2  *C0*C4)  * M3* (01**2-2*00*02)  *N4+C0**2*M5)  / (2*DEL6)  i 

go  to  10  V 'r:-Sn 

7 CONTIN  02  " r - , • ! 

C*********** *************************** ********************************< 

c**  *< 

C**  THIS  SEGMENT  CALCULATES  THE  SEVENTH  ORDER  SYSTEM  INTEGRAL.  **. 

C**  ' • V ---  > - *< 

£*****  ******  ******  * ***********  * *****  ********  *******  ***********  *********  4 

READ(5, 107)06,05, C4, C3, C2,C1,C0  _ : ‘.i.  i-:| 

READ  (5,  108)D7,D6,D5,  D4, D3,D2,D1 , DO 

N1=-(D1*D4-D0*D5)  **  2*  (D0*D3-Dl  *D2)  *(  DO  *D7-D  1*D6*D2*D5-D3*D4) 

(12=  (D0*D7-D1  *D6)  * (-DO  *D5+D1  *D4)  ♦ (D0*D3-D  1*D2)  * (D2*D7-D3*  D6)  \ 

H3=-  (D0*D7-D1  *D6)  **2+  (DO  *D3-D1  +D2)  * (D4*D7-D5*D6). 

(14=  ( 1/DO)  *(D2*M3-D4*M2+D6*M1)  ':r‘v;€W- 

M5=  (1/DO)  * (D2*M4-D4*M3+D6*M2)  . /v.: ' . ±- 

(16=  (1/DO)  *(D2*M5-D4*M4+D6*M3) 

M0=  (1/D7)  * (D5*M1-D3*M2*D  1*M3) 

DEL7=D0*  (D1*M6-D3*M5  + D5*M4-D7*M3)  .......  _ i..._ 

BNTGRL=  (1/(2  *DEL7) ) * (C6**2*M0*  (05**2-2*04+06)  *Hl+  (04**2- 2*C5+C3*2* 

; 102*06)  *M2*(C3+*2-2*C2*C4*2*C1*C5-2*C0*C6)  *M3 ♦ (C2**2-2*C  1*C3+2*C0*C 
--24)  *H4+ (C  1**2-2*00*02)  *K5+C0**2*M6)  ■ 

GO  TO  10 

8 CONTINUE 

C*** ********  ******  ****************************************************** 

c**  - - ■ -5~- -¥§£■ >?•' 

C**  THIS  SEGMENT  CALCULATES  THE  INTEGRAL  FOR  AN  EIGHTH  ORDER  SYSTEM** 

C**  ' ..  . ..  ** 

C** *********************************************************** ********** 

READ  (5, 103)  07,06,05,0  4,0  3,02,01.00 

READ  (5,105)  D8,  D7,D6,D5,D4  ...  

READ  (5,  104)  D3,D2,D1,D0  ■ 

Ml=(D0*D7+D2*D5)  * (-D0*D1  *D7+ D0*D3*D5+2 +D 1 **2*D6)  ♦ ( D3*D7-D5**2)  * 

1 (D0**2*D5+Dl*D2**2)  +D  1*D 3*D 3*  (D 0*D3-D1  *D2)  -D1**2*D3*  (D0*D5-D1  *D4 ) * 

2 (-D2*D7+D3*D6-D4*D5)  * (DO  *D3  **2  ♦ D1  **2*D4)  -D1*D6*  (D1**2*D6* 
33*D0*D3*D5)  -D1  *D2*D3*  (D3  *D6-D4*  D5)  +2*D0*  0 1 + D4*D5+*2 

M2=(D0*D3-D1*D2)  ♦ (D0*D7**2-D1*D5*D8-D1*D6*D7^D2*D5*D7)  + _ I 

1 (D3*D8-D4*D7)  * (-D0*D1*D5  + D0*D3**2-D1*D2*D3*D1**2*D4)  - : % . 

2D0*D5*D7*(D0*D5-D1*D4) +D1**2*D8*  (D0+D7-D1 +D6)  : 


1-D1**2*D4)  +P0*D7**2*  (-DO  *D5 + D1  *D4+D2*D3)  - (2*D0*D1*D3*D7*D8) 

(14=  (-D5*D8  + D6*D7)  * (2*D0*D1*D7-D0+D3*D5+d1*D2*D5-D1**?*D6)  ♦ 

1 (-D3*D8+D4*D7)  *(D0*03*D7-Dl*D2*D7+Dl**2*D8)-D0**2*D7**3 
H5=(D2*(l4-D4*M3  + n6*M2-D8*Ml) /DO  . : :Xv 

(16=  (D2*M5-D4*il4*n6*M3-D8*M2) /DO 

(17=  (D2*M6-D4  *M5+D6*M4-D8*M3)  /DO  . :^V?v  'Jl. 

(10=  (D6*M1-D4*M2+D2*M3-D0*M4)  /D8 
DEL8=D0*  (Dl*tl7-D3*M6  + D5*  M5-D7*(14) 

ENTGRL=  (C7**  2*M0*  (C6**2-  2*C5*C7)  *M1+  (C5**2-2*C4*C6>2*C3*C7)  *H2* 

1 (04**2-2*03*05 ♦2*02*06-2*01 *07) *M3+ (03**2- 2*C2*C4*2*C1*C5-2*C0*C6) 
2*H4*  (C 2* *2- 2*0 1*03+ 2* CO *04)  *M5*  (C  1**2-?  2*00*02)  *M6*C0**2*H7) /(2+DEL 

38)  - • ...  - ■- 

GO  TO  10 
9 CONTINUE 

Q* **********  ******  ****************************************************** 

c**  " ** 

C**  THIS  SEGMENT  CALCULATES  THE  NINTH  ORDER  SYSTEM  INTEGRAL.  *♦ 

C**  +* 

£**«** *****************  *******************  ***************  *************** 

READ (5, 105)  08,07,06,05,04 

R3AD(5,  104)C3,C2,Cl,C0  ...  ^ .. 

ice  ‘ 


READ  (5,1 05)  D9,  D8,  D7,  D6,D5  v-v.-vif-fe'.  , • ,=4i  v : f 

READ  (5,105)  D4,03,D2,D1,D0  f ^ 

Al=D1*D2-D0*D3  -----  . •-  - - j| 

A2=D1*D4-D0*D5 

A3=D3*D4-D2*D5 

A4=D1*D6-D0*D7  r -vs- — - --r  I 

A5=D3*D6-D2*D7  . j ,y  - - ^1;!^  '•  ' > l 

A6=D5*D6-D4*D7  : - •;•; ;,v  >v  :,v  J 

A7=D1*D8-D0*D9  j 

A8-D3*  D8-D2  * D9 
A9=D5*D8-D4*D9 

A10=D7*D8-D6*D9  . " ...-.•  - . , . ■ rv.^  r ■:..T.;rv:-a«r.  -s«C.,3rSi8a,vsS.! 

y-  Hl=A1*  (A1*A1  0-A2*A9+A3*A6+A3*A3+2*A4*A6-A5**2-A5*A7-;i7**2)  £ 

1A2*  (-A2*A6-A3*A7+A4*A5+2*A4*A7) -A4**3  V-:  T W-'-  * ' 

H2=A1*  (A3*A9  +A4*A9-A5*A8  + A5*A7-A7*A9)  + A2*  (-A2*A9*A4  *AB+A7**2)  -A4** 
12*A7 

H3=A1*  (A3*A10+A4*A10+A7*A9-A8**2)+A2*(-A2*A10+A7*A3)  -A4*A7**2 
- H4=A1*  (A5*A10+A7*2*A10-A8*A9)  + A2*(A7*A9-A4*A10) -A7**3  ‘ 

85=  (D2*M4-D4  *8  3*'D6*M2-D8  *H1)  /DO  • - . • 

. H6=(D2*W5-D4*fI4  + D6*M3-D8*H2)  /DO  " y'r  - - ty 

H7=  (02*86-04  *M5+D6*H4-D8  *83)  /DO  -JiSeBUL-sSasSU 

88=  (D2*M7-D4*H6+D6*M5-D8*B4) /DO 
H0=  (1/D9)  * (D7*M1-D5*H2+DJ*83-D1*34) 

. DEL9=D0* (D1 *H8-D3*H7 +D5*H6-D7*M5+D9*H4 ) , -rr-,  T- 

ENTGRL=(C8**2*fl0  + (C7**2-2*C6*C8)  *81+ (C6**2-2*C5*C7+2*C4*C8) *82+  , j 
1 (C5**2-2*C4*C6+2*C3*C7-2*C2*C8)  *83+(C4**2-2*C3*C5+2*C2*C6- 
22*C1*C7+2*C0*C8)  *84+ (C3**2-2*C2*C4+2*C1*C5-2*C0*C6)  *85  . \ 

3+  (C2**2-2*Cl*C3+2*C0*C4)  *86+  (Cl  **2-2*C 0*c2)  *B7+CC**2*M3)  /(I*DBL9) 
C***** *********************************************** *******************’ 


>85  10  CONTINOE 

86  ENTGRL=ENTGRL*  2*  3. 1 41  59 

87  ENTGRL=ENrGRL**0.5 

88  WRITE(6,301)  (KE30  (MIN)  ,NTNS 

89  WRITE (6, 200)  EHTGRL 

90  999  CONTINUE 

91  STOP 

92  END 

SENTRY 

THE  BBS  VALUE  IS  GIVEN  FOR  EACH  INPUT. 


1,15) 


**  NOTE: 

0.32670-02 
**  NOTE: 

0.2650D  00 
♦*  NOTE: 

0.  178  OD  00 
**  NOTE: 

0. 18150-01 
••  NOTE: 
__0.8 19  00-03 
MOTE: 

0.21240-02 
••  MOTE: 

, 0.32660  00 
T>*  NOTE: 
>0.21780  00 
NOTE: 

0.24990-01 
**  NOTE: 
0.550SO-03 


OPTIBAL  C.3.  ACCELERATION,  RUN  B4 
TRAILING  EDGE  RHS  VARIATION 


LEADING  EDGE  RHS  VARIATION 


OW1  RHS  DEFLECTION 
DW2  RHS  DEFLECTION 

OPTIBAL-  C.  G^VCCELER ATIO Ny  ITOM-? B5-:=g 
TRAILING  FOGS  GAP  RHS  VARIATION 
LEADING  EDGE  GAP  RHS  VARIATION 


DW1  RHS  OEFL ECTION 


ON 2 BHS  DEFLECTION 


v dtx  pic  l 

(1J  * OPTimL  PROBLEtU  . * liPl ) * . * C«N5> 

(21  N-  1**41 
(31 

(41  Z2«-01«**l 

(53  r«*zt,mz2 

(6)  theta-dti  sm-z 

(73  THET*U+lN,m*THETP  “ 

(03  TH£TPl2~iN.-NitTHETP 
(93  THEr*2l»H-N),N)rTHET* 

( 10  3 THET&22-K  -N,N  it  THETA 
( 113  R+iN*N)>  0 
(123  1-0 

C 13 1 LOOP  l POL  D+P 

(14  3 P-t  THETH2 1 * THE  TP22  ♦ . */?>  ♦ .«•<  THETP11*TH£TP12*  . */?> 

(153  P*<ff*«*>«2 

(163  EPPQR~<.*/*/lROLD'mR)  *2)  *0*5 

1171  I-I*l 

(103  CPPOP 

(191  MERPOR>D 'LOOP 

(201  /cr»l0Pl>*.»<(«M5)*(<«01)*.*ff» 

(211  ‘Pr-* 

(22)  io  "3  wr 

(231  R+RU-BU*.*KT 


104 


I 


111 

121 

131 

14] 

15) 
161 
t7) 
LSI 
C 9 ) 

r ioj 
mi 
I 12) 
I 13) 
I 14) 

(15) 

(16) 


Cl) 

C2) 

(3) 

(4) 
(«) 
16) 
17) 
C0) 
19) 
CIO) 
Cll) 
C 12) 
(13) 
(1^1 


v COSTFN 
Z-6*0 

•ENTER  Cl  THRU  C6' 

z-a 

0+  6 6 '0 

_ _ _ 

0-  0 1 234340 

6 7 »0 

tt+TR ... 

05«-t7,7»»  ALL 

H5+l?,-2>*ALL 

Pl-t-2,2)tfiLL 

01- 05- CN5* .* CQPl >♦.* CHN5) ) 
•you  Mfff  MOM  READY  TO  RUN  RIC% 
Al+AU-BU*  . * CiPl ) ♦ . » C«M5> 

Bl+BU 

9 

WJMSCO)* 

V UlNGlZlIlP 


C15)  P+POt  14 ) 

CIS)  LOOPEH+l'l 

(17)  ••(UPCIllilZ  6)SL00P2 

CIO)  X*C»P)*W 

Clf)  P^C-DtP 

C20)  fP 

7 


* THE  DUX  POLYNOMIAL  ISt • 
J*0 

P*P0(  )3) 

LOOPtl+I*l  _ 

Mi  \PinilUT6) SLOOP 
IMW1-Z 
P»<-I)tP 
P P 

1 THE  DUE  POLYNOMIAL  IS* 1 

XM 


▼ SEAT  XO  iPAlPBlAUAZ 
(11  41-4 

(21  Pit  *11— * 

(31  42-41 

(41  42(1*  1*0 

(81  PA+iFRAME  AD-xFRAME  *2) 

(61  41-4 

(71  41(*21— B 

(81  42-41 

(91  42(2*1-0 
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V GAPiC0iPlliJ;K*Al}A2 

I 1 1 ft  THIS  routine  calculates  the 

12  i •L.E.  AND  r.£.  CAP  NUMERATORS 

13  J •AND  THE  RINGLET  MOTION 

14)  * numerators . 

(5)  1 f. these  are  to  be  used  to  find 

(6)  •the  RESPECT IOE  RMS  VALUES. 

17 J * THIS  FUNCTION  CALCULATES  THE 
183  ^LEADING  AND  TRAILING  EDGE  GAP 
193  •NUMERATOR  POLYNOMIALS . 
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•make  l u l?  til 

v rtAkEtNtff  (21 

ill  ' A GENERATES  PASSIVE  a matrix  (31 
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V 0PT}J0LDtZ2tI tZSTEPiN 
(11  *TH1S  PROCEDURE  PERFQRHS  A 
(£1  *FOUR  UARIABLE  GRADIENT  SEARCH 
t 3 1 UNROUND  ft  FUNCTION  OF  COST 
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